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ABSTRACT 


Problems of control are presented from a standpoint of simple 
operation rather than mere mechanical simplicity. Of first im- 
portance is the elimination of inconsistent or abnormal behavior 
in which spinning is classed. The influence of wing and engine 
position is discussed, and several advantages of an automobile- 
type wheel arrangement are brought out. 

A primary control, in which flaps are coordinated with elevator 
and rudder with ailerons, is advocated. In a possible supple- 
mentary pedal control, ground brakes are coordinated with air 
drag obtained by opposite deflection of two rudder surfaces. 

Special problems, considered from a standpoint of simplified 
control, include stability against overturning on the ground, 
inertia loads on the wing, spiral stability, and control require- 
ments for following a compass course. 

The ideal airplane for private flying is described as being: out- 
standing in vision, incapable of spinning, comparable with an 
automobile in simplicity of control, yet with unquestioned 
superiority of cross-country performance. 





oo CONTROL, from the pilot’s standpoint, must 
have its foundation in positive and consistent 
relations between the control application and its re- 


sult. This requirement bans any abnormal flight 
condition which, even for a brief interval, amounts to 
loss of control—like the stalled landing still so com- 
monly in use. The development of a more consistent 
landing procedure, although started where most im- 
portant, with private planes, has now been applied to 
almost all types. Other improvements will not long be 
denied. For example, instead of studying spins as 
normal or abnormal, the broader and sounder view is 
now being adopted—that any spin is abnormal and 
hence should be outlawed from the physically possible 
maneuvers of any airplane designed for transportation. 

Although the nonspinning airplane has only lately 
been officially recognized, it is not particularly new or 
difficult of attainment. The earliest such plane known 
to the author which was essentially in this class was the 
Aeromarine-Klemm.! The original design as imported 
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from Germany was so far from being nonspinning that 
it would, on rather slight provocation, go into a flat 
spin from which recovery proved impossible. For- 
tunately, the light wing loading minimized the fatalities 
from several such occurrences, but still it was no plane 
to put in the hands of a novice. ° Hence, the aero- 
dynamic properties were thoroughly revised, the modi- 
fications including less taper, a thicker tip section, less 
control movement and further forward c.g. In the 
absence at that time of precise knowledge as to the 
effect of such changes, most of them could be classed as 
mere engineering judgment, and at least one—the re- 
duction of taper—was quite fortuitous. The result was 
most gratifying. After test pilots were unable to spin 
the plane, the air commerce inspector finally succeeded, 
by a combination of whip stall and crossing of controls, 
to hold about two turns, after which it would recover 
by itself. This result seemed even better than being 
totally unable to spin, because the tradition was then 
quite strong that when a ship was hard to spin it was 
also difficult to pull out. 

The next step in the development of simplified control 
may well be ascribed to the first Stout Sky Car of 
eleven years ago.?, Again there was no attempt to make 
it absolutely nonspinning but merely spin resistant to 
the extent that it could not fall into a spin inadver- 
tently from any conceivable attitude of normal flight. 
It had essentially the modern tricycle landing gear, 
except for the steerable front wheel, and was also de- 
signed for what is now known as two-control operation 
(although the rudder pedals were retained for those 
who insisted on using them). As originally assembled 
with the floating wing-tip ailerons, it had the most per- 
fect lateral control of any plane the author had ever 
flown. By overbalancing the ailerons with weight for- 
ward they were put normally in a position of slight 
negative lift, which in turn gave a slight yaw in the 
favorable sense when operated. As a result, almost any 
degree of turn could easily be made with a quick entry 
and a perfect bank by using the stick alone and with 
feet entirely off the pedals. 
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Either of the above-mentioned aircraft, and doubt- 
less others as well, would have been completely spin- 
proof if independent rudder control had been elimi- 
nated. 

Further developments along similar lines have in- 
cluded the Weick W-1, the Hammond Y, the Gwinn 
Aircar, the Ercoupe and the General Skyfarer, to name 
only those airplanes in the unconventional class de- 
signed for private use. Most of these have shown an 
unmistakable trend toward spin-proof qualities, the 
level-landing type of wheel arrangement, front wheel 
steerable by the control wheel, two-control operation 
and a somewhat hesitant use of flaps. Questions that 
are still unsettled are: pusher versus tractor, high-wing 
versus low-wing, important features of the landing 
gear, and specific means of operating and coordinating 
the various controls. 


ENGINE AND WING POSITION 


Although the position of the engine and wing is 
usually dictated by considerations other than control, 
the following control factors should be kept in mind: 

The pusher engine involves more shift of center of 
gravity as between full load and light load, but, on the 
other hand, can be arranged to give less change of trim 
than a tractor as between power on and power off. One 
reason is that the pusher engine can be mounted con- 
siderably higher, thereby helping to counteract the 
usual nosing up tendency when power is applied. An- 
other is that more effect can be obtained by a change in 
slope of the thrust line in a vertical plane, this being of 
relatively minor effect in a tractor because of the pre- 
dominating influence of the wing on the direction of the 
slip-stream. 

Maintenance of the same stalling characteristics for 
power on or off is also a distinct contribution of the 
pusher engine to simple, consistent control. Although 
the slipstream effect on the wing lift is thereby lost, the 
time when such help would be most welcome is usually 
when it is nonexistent. 

The main advantage of the pusher .is, of course, 
vision, the all-important reference on which any control 
must normally depend. The pusher arrangement is 
particularly desirable if a high wing is used. 

As is well known, the high-wing position makes both 
longitudinal and lateral stability easy to obtain, es- 
pecially at the more critical high angles of attack, lead- 
ing to a type of control that is generally more positive 
and consistent. This does not mean that simple and 
satisfactory control cannot be had with a low-wing ship 
but only that such a case is the exception rather than 
the rule. 


TRICYCLE LANDING GEAR 


As the landing and, to a lesser extent, take-off have 
hitherto been the most critical parts of a pilot’s primary 
education, the landing gear calls for first attention in 
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the discussion of control means. Without going into 
the now well-recognized advantages of the so-called 
tricycle gear,’ certain faults in its operation, and whether 
there is not still room for considerable improvement, 
should be considered. 

One criticism that has been seriously made of the 
tricycle gear is this: Any pilot placed in such a plane, 
even though experienced and generally sensible, after a 
period of somewhat timidly feeling out its character- 
istics and satisfying himself that it can be landed at the 
will of the pilot instead of at the will of the plane, will 
adopt the attitude that he no longer needs to be careful. 
The plane will be landed with little regard for the wind, 
with practically any speed of descent, and using the 
brakes without mercy until something gives way. This 
indicates the need for unusual attention to strength 
factors and brake capacity if the inherent advantage 
made mechanically possible by the tricycle gear is to be 
largely realized. With respect to side-wind landings, 
there is, at the same time, a rather basic factor that has 
apparently not yet been given adequate attention. 

A primary condition for simple consistent control is 
obviously that the control movement in a given direc- 
tion must always serve a similar function. Examples 
of inconsistent control, familiar to any pilot of a con- 
ventional plane, include pushing forward on the control 
when the nose goes down in a spin, holding the control 
back at the moment of landing, and the formerly com- 
mon practice of moving the aileron control away from 
the side toward which it is desired to turn on the ground. 
The first example is eliminated in the nonspinning 
plane; the second, by the tricycle or level landing gear; 
while it might be supposed that the third would be 
eliminated by the steerable front wheel geared to the 
control wheel. Unfortunately, however, conditions are 
thereby introduced which may make the aileron con- 
trol still less consistent. 


STEERABLE FRONT WHEEL 


Whether two-control or three-control operation is 
used on the immediate approach to a landing, the wings 
must be brought to a reasonable level by use of the 
control wheel. Suppose that as a result of a cross-wind 
landing or offset position of the control, one wing per- 
sists in going up after ground contact has been made. 
To be entirely consistent, the pilot would then attempt 
to lower the high wing by rotation of the control wheel 
in that direction. This, however, turns the front wheel 
now in contact with the ground, in a direction to throw 
the centrifugal force the other way, which invariably 
overbalances the aerodynamic effect of the ailerons and, 
instead of correcting the rolling tendency, actually 
makes it worse. 

Consider the example of a pilot approaching for a 
landing and discovering that he is drifting laterally as 
he nears the ground. If there is still room for a banked 
turn, he, of course, uses the remaining height to turn as 
far as he can into the wind. On closer approach to the 
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ground, however, he brings the wings level, and to stop 
them in the level condition he reverses again momen- 
tarily. In this way actual contact may be made with 
some of the same lateral drift still remaining and with 
the front wheel turned somewhat into the wind. This 
puts the centrifugal force, as well as the primary lateral 
drift force, in a direction to cause a down-wind over- 
turning tendency, with the resultant inconsistency of 
control reaction already mentioned. If a wing is 
damaged it may be called “‘pilot’s error,” but the fact 
remains that there are sure to be plenty of such errors 
in any airplane that invites them. 

A fixed front wheel seems definitely better than one 
that might be turned wrong. Yet with the steerable 
front wheel accepted as generally desirable, there are at 
least three cures for this tendency to steer in the wrong 
direction: First, let the pilot forget that he is flying an 
airplane as soon as he makes contact with the ground, 
and then imagine that he is in an automobile to be con- 
trolled accordingly (by steering into the bank). Sec- 
ond, with a suitably designed front wheel, simply re- 
lease the control entirely and: let the dynamics of the 
situation take charge.* Third, provide a wheel spacing 
that in itself will prevent overturning. Neither of the 
first two measures is likely to be sufficient in itself be- 
cause both depend on abnormal action or inaction by 
the pilot. Even aside from wrong control, the possi- 
bility must be considered of the front wheel dropping 
into a rut or that positive steering must be used to 
avoid obstructions. Hence it is important to provide 
enough wheel base and tread to take care of any situ- 
ation that may be reasonably expected. This condi- 
tion is particularly critical with the high-wing arrange- 
ment because of its relatively high center of gravity 
and also because the most logical mounting for the 
wheels in that case is on the fuselage or body. 














a’ 
1 : 
od 


CH Empty 





Fic. 1. Landing force diagram. 


DEVELOPMENTS IN SIMPLIFIED CONTROL 


THREE VERSUS FouR WHEELS 


Fig. 1 shows two views of an airplane body in which 
the dotted wheels show a tricycle arrangement and the 
full-lined wheels show the equivalent four-wheel gear. 
The arrow in each case represents the weight and in- 
ertia force projected into the plane of the paper, with 
the brake-force coefficient 0.55 and the lateral force 
0.75 of the vertical component. The solid arrow (0.93 
resultant ground force) represents the condition with 
brakes applied on all wheels; the dotted arrow, with 
brakes on the rear wheels only. With the scale chosen 
(which makes the force vector just reach the ground), 
the limiting condition at overturning is represented by 
that which makes the point of the arrow in the top 
view just touch the lateral wheel line. It is evident 
that with all wheels brake-equipped the rear wheels of a 
tricycle gear must be placed at an almost prohibitive 
distance from the body to get the same ground stability 
against overturning as is provided by the four wheels 
shown, the latter having no greater tread width than 
that of an automobile. The total brake force obtain- 


able from all wheels is over 60 per cent greater than 
with brakes on the rear wheels alone, the relative stop- 
ping distance being approximately in inverse propor- 
tion. While there is some question as to the validity of 
starting with the two ground-force components rather 
than with the magnitude of the force itself and its 
longitudinal component, the conclusions are qualita- 


tively the same in either case. The term “ground 
force’’ refers to any force in a horizontal plane trans- 
mitted by the airplane to the ground. 

Contrary to automobile practice, it seems desirable 
to give the front wheels a tendency to lock slightly be- 
fore the rear wheels, although the locking of any wheels 
should be possible only with a very substantial pedal 
force. In general, the ground stabilizing factors are 
chosen with a view of resisting a turn rather than 
accentuating it and without much concern for tire wear, 
which is normally of small importance in an airplane. 
It may be remarked in passing that the commonly used 
circumferential ribs on airplane tires are quite unsuited 
to either type of landing gear. Whereas an important 
object in an automobile is to prevent skidding, the ob- 
ject in an airplane should be to permit ground skidding, 
particularly of the front wheels, at any lateral force 
approaching a critical value. 

Because of the offset position of the kingpin pivot 
axis for each of the two front wheels of the four-wheel 
gear, there are other important geometric variables in 
the linkage besides the caster angle utilized in the single 
front wheel. For example, proportions can be readily 
chosen such that the stabilizing tendency of the wheels 
will be accentuated by brake application or soft ground. 
For means of avoiding ‘‘shimmy,” there is a wealth of 
experience in automobile practice. In a high-wing air- 
plane of a size permitting a choice of suitable wheels, 
there is little drag or weight difference between the 














518 JOURNAL OF THE AERONAUTICAL SCIENCES—DECEMBER, 1942 


three and four-wheel gear; the latter affords an appre- 


ciable margin of security against the failure of a wheel 


or tire, particularly in front. 

Advantages and special characteristics of the four- 
wheel gear have been brought out not by way of arguing 
its adoption in all cases but to show what can be done 
with an alternative arrangement. In actual practice, 
the number of wheels, like the position of the wing, 
must depend largely on considerations of purpose and 
general arrangement. 


POSITION OF REAR WHEELS 


If the tail is to be designed for no ground contact and 
if the plane in general is to resist a reasonable ground 
wind without staking down, it is important to put the 
rear wheels a substantial distance behind the most rear- 
ward c.g. To a certain extent this requirement fits in 
with the desire to prevent taking off prematurely. Yet 
there must be no interference with getting a sufficient 
angle of attack for a proper take-off at a speed prefer- 
ably about 20 per cent above stalling. The combined 
requirements are most difficult of attainment for a 
high-wing pusher because of the high center of drag and 
weight and the large rearward shift of the c.g. on 
removal of disposable load, not to mention the thrust 
force. The original Sky Car, for example, would im- 
mediately tip back on its tail when the occupants got 
out. This action was taken care of by having a tail 
wheel in addition to the nose wheel, but at the best it 
was disconcerting to a passenger. At the worst it made 
the ground attitude of the empty ship no better than 
with a conventional landing gear as far as reaction to 
ground winds was concerned. 

If at the moment of take-off the effective angle of 
attack can be increased sufficiently by flaps, it is 
obvious that a take-off can be made with the rear 
wheels in almost any position. But, although flaps are 
useful in take-off for other reasons as well, the use of a 
separate flap control for the purpose would be anything 
but a simplification. The solution of this problem can 
be combined with other desired results in the method 
of coordinating the longitudinal control. 


Lirt CONTROL 


The primary lift control has always been, and will 
continue to be, the fore-and-aft movement of the con- 
trol column, traditionally through the medium of the 
elevator. The only other instrumentality involving 
lift in any major sense is the flap. Hitherto, the flap 
has usually been operated as a combined lift and drag 
control by means of a separate flap control lever. Such 
an arrangement is admittedly simple mechanically 
but complicated in operation, since operating simplicity 
requires coordination of like rather than unlike func- 
tions. Hence, with a relatively low-drag flap’ (of which 
there are several types now available), the rational 








place for its control is in direct coordination with the 
elevator. 

In order to minimize for the normal cruising range 
such drag as is still unavoidably involved in flap opera- 
tion, the flap movement should be designed to pick up 
in parabolic relationship to the elevator as the control 
column is pulled back. The effect on the wing section 
characteristics of the combined control in comparison 
with the normal elevator is shown in Fig. 2. In the 
lower part of this figure the heavy full lines show the 
increase of lift coefficient with angle of attack of the 
basic airfoil (and of the airplane bodily) as produced by 
operation of the elevator alone. The heavy dashed 
line cutting across the lift curves for varying flap de- 
flection shows the effect of the proposed flap and 
elevator combination. The corresponding drag is 
shown above. 
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Fic. 2. Example of coordinated lift with incidental"drag. 


With a proper linkage ratio in the system, there 
appears no reason to expect an essential operating diff- 
erence in the process of getting lift whether by angular 
deflection of the entire airplane or by combining this 
process with that of deflecting the wing or a part of it. 
The detail characteristics of the coordinated control 
are subject, however, to a few qualifications that include 
the following: 

(1) The direct use of the flap facilitates take-off not 
only by reason of the lift itself (which can be withheld 
until the moment when it is needed) but by permitting 
the application of substantial lift without depressing 
the tail, as already mentioned. 
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(2) Clearing an obstacle is made somewhat more 
positive by the larger lift quickly available, as well as 
by the fact that for a given lift the tail need not be 
depressed so much. 

(3) For given balance characteristics, the control 
forces are increased by an increment that is desirable 
in the smaller planes, but in much larger sizes and higher 
speeds will undoubtedly set a somewhat lower limit to 
what can be handled without boosters or servos. 

(4) There will be more aerodynamic margin of safe 
control in maneuvers such as sharp turns and pull-ups 
because, for a given air attitude of the airplane, a 
greater effective camber of the wing is automatically 
provided to carry the load put on it. 

(5) The control reaction will be so direct and un- 
escapable that there would appear to be little or no ex- 
cuse for inadvertent mistakes such as too suddenly 
retracting the flaps. It must be expected, however, that 
the airplane lift will always react more quickly to the 
flap than to elevator deflection. 

(6) For the wing stress analysis, instead of con- 
sidering flaps-down as a supplémentary condition, each 
major condition is governed by its own flap deflection. 

(7) A maneuvering load factor can be realized more 
quickly but with a more substantial control force, so 
that on the whole there would appear less reason for 
inadvertent overloading. 

(8) Landing loads should be appreciably less 
because of the facility of quickly checking the sinking 
speed without waiting for the tail to drop. Like any 
other plane, however, the glide-in should never be made 
at the rearward control limit. 

The feasibility of gliding in with a substantial speed 
reserve, at the same time maintaining a steep approach, 
depends, of course, on drag as well as lift. 


DracG CONTROL 


There is no longer much room for argument on the 
desirability of separating drag as far as possible from 
lift and employing it as a secondary control.‘ The point 
is well realized by any pilot who has dropped down 
toward a small field with split flaps set at a large angle 
only to find, on close approach, that the field is unusable 
and that he must climb out again. 

With the object of utilizing other existing surfaces 
and controls if possible, the following appear to be the 
most plausible means of effecting the desired functional 
separation, alternatively: 

(1) Design the flap so that beyond the lift range 
there is a range of sharply varying drag but practically 
constant lift. This part of the movement obviously 
should not be kept on the main control column but 
might be put on a separate lever. This still, however, 
leaves several apparent chances for confusion of control. 

(2) Use two rudders or a split single rudder so that 
the two surfaces may be deflected outward by pressure 
on a foot pedal. There must be at least one pedal for 
brake operation on the ground, and it can be used for 
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the aerodynamic brake as well. The use of the rudders 
provides ample drag for glide path control (up to a 
gliding angle of about 1:7); they can be readily ar- 
ranged to retain a differential movement for lateral 
control purposes and, by their position at the rear, give 
an added stabilizing effect both directionally and longi- 
tudinally. 

In a broad sense, directional control can be con- 
sidered as differential drag; hence, any ship for which 
some measure of independent directional control is 
desired can be equipped with the conventional two 
pedals instead of one and can be arranged so that the 
depression of one or the other pedal turns one or both 
rudders in the corresponding direction. The application 
of a symmetric drag force will then be by depressing 
both pedals. Somewhat similar independent action of 
the two pedals has been used before (notably on tailless 
airplanes) and was tried out by the author several 
years ago with complete success. The fear expressed 
by some that there might be danger of inadvertently 
pressing on both pedals at an inappropriate time, such 
as during the take-off, did not seem to be borne out in 
practice. In fact one flight instructor thought the 
possibility of such an occurrence was a distinct training 
asset because it would so positively discourage undue 
pressure and mental tension on the part of a student. 

For most purposes, of private flying at least, it 
appears preferable to stick to the single pedal and 
coordinate the usual rudder control entirely with the 
ailerons. In the interest of functional separation, it is 
important to position the rudder surfaces in such a way 
that their deflection will have no pitching effect. 


STEERING CONTROL 


Coordinated aileron and steering control is, of course, 
not essentially new. The original Wright airplane 
employed it in principle, and almost any well-designed 
modern plane can be steered in a reasonably safe man- 
ner on the ailerons alone. There appear to be good 
reasons, however, for designing a definite two-control 
airplane on the basis of direct coordination of actual 
rudders with the aileron system (as did the Wrights) 
rather than by relying on sideslip for supplying the 
necessary lateral tail.reaction. 

Two-control ships hitherto available have called 
forth considerable criticism from trained pilots who did 
not feel at home in a plane that yawed and slipped per- 
ceptibly and yet which gave them no means by which 
it could be remedied. Such opinions, based on flight 
habits in a conventional type of airplane, mean nothing 
in themselves—at least from the standpoint of safe 
control. The danger of a slip, and particularly of a skid, 
has nothing directly to do with the maneuver itself but 
only with the dangerous attitude of the airplane that 
such a maneuver may lead to; any plane that is pre- 
vented by other means from staying out of a spin is still 
safe no matter how much it may slip or skid in the turn. 
Such criticism as the above has usually been applied to 
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relatively sharp turns where it is really not warranted. 
Experience, however, has indicated a real objection to 
departures from a perfect bank as applied to slight 
turns. 

This objection has to do with the procedure of follow- 
ing a compass course for which the flight path, instead 
of being straight, is in reality a series of flat curves. A 
plane that depends on slipping for its directional control 
takes an appreciable time to swing into the turn after 
banking, and it similarly skids and wants to keep turn- 
ing on the recovery. Hence, the tendency is for the 
plane to be late going in and late again coming out so 
that it turns beyond the heading anticipated. It hap- 
pens that the compass on a southern heading acts just 
opposite to this tendency and encourages straightening 
out too soon; hence, on this heading, the control and 
instrumental tendencies supplement each other, with a 
result that is usually satisfactory. On a northern head- 
ing, on the other hand, the two tendencies add up to a 
result that is worse than either and which by actual ex- 
perience with an ordinary type of compass can be so 
serious as to make difficult the following of a reasonably 
straight course in even mildly rough air. Experience 
shows, however, that it is readily possible to make a 
rudderless airplane behave better than one with im- 
properly coordinated rudders.§ 

As shown by Fig. 3, even differential ailerons give un- 
favorable yaw moments at small deflections.6 To 
remedy this condition, the coordinated rudder move- 
ment should be enough to overbalance substantially the 


Gs pA 





Rudder|CG"x/0 





















































7 
v ws en waged 
Q 0 me 20° 
% VA i Aile Ci x10 
Ss a 
g a ; 
rs eral 
w J] at : 7- 
2 rs be 
+~ 
a ee ee 5 (in) — 
A OE Ys 
as 146 jar / 20. Bere" 
Alleron Deflection 
a (degrees. 
Q 
702 
S eA 
& —— 
S (Moment coerricients from Ref.6) 
| -5*Up aileron deflection 
& *Down “ cs 
5, =Rudder # (outward & inward) 


Fic. 3. 


Example of coordinated steering by two ailerons and 
two rudders. 
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aileron yaw, particularly at small deflections, the same 
as would be done consciously or subconsciously by a 
good pilot under similar circumstances. The proper 
magnitude of directional overbalance depends, of 
course, on the relative rolling velocity and the moment 
of inertia around the vertical axis. Another considera- 
tion is the unfavorable rolling arm of the usual rudder, 
which should be kept as small as possible. 


SPIRAL STABILITY 


A disadvantage of the rudderless airplane is the spiral 
instability necessarily involved. It is true that prob- 
ably no airplane flying today is spirally stable to any 
appreciable or practical extent, regardless of its control 
system. This condition is not so much inherent, how- 
ever, as due to the somewhat arbitrary restrictions of 
dihedral and other characteristics imposed by the con- 
ventional control system. For example, aileron control 
by itself quickly becomes indeterminate, or even has a 
reverse effect on roll, as the dihedral is increased. With 
adequate coordinated rudder movement, however, sub- 
stantial dihedral will help rather than hinder the attain- 
ment of the proper bank. 

The amount of dihedral necessary for positive spiral 
stability can be minimized by resort to other expedients 
such as the introduction of effective fin area forward to 
supplement a corresponding increase at the tail in re- 
sisting yaw and sideslip. One of the best ways of doing 
this is by use of a compact streamlined body of rounded 
section such as will give generally efficient performance 
in a pusher type of plane.’ Such a body has its lateral 
center of pressure far forward so that it is equivalent to 
a fin surface considerably forward of the nose end. An- 
other possibility of promoting spiral stability is by 
direct damping of turn by drag elements near the wing 
tips which will automatically function by differential 
tip speed. 

Any model builder knows that it is possible to achieve 
spiral stability with nothing more than adequate di- 
hedral, although this by itself will probably not be an 
entirely satisfactory means for a practical airplane 
(largely due to the effect of side gusts).® Still, with so 
many possibilities available the author does not fall in 
with the commonly accepted attitude that substantial 
spiral stability is an impractical dream. The solution 
of this problem by simple means would at one stroke 
remove the major difficulty and danger from blind 
flying. A good illustration of this fact is that airships 
of the lighter-than-air type have never imposed any 
problem of control in blind flying. The obvious reason 
is that they are spirally stable, even though less stable 
than most airplanes in a purely directional sense. 


CONCLUSIONS 


Although it is believed that many of the features de- 
scribed can be adopted advantageously for other types 
(certainly including the military), they are most di- 

(Continued on page 548) 








Isaac Newton and Aerodynamics 


TH. VON KARMAN* 
California Institute of Technology 


INTRODUCTION 


O* CHRISTMAS DAY of the year 1642 in a manor house 
in Woolsthorpe by Colsterworth in Lincolnshire, 
England, a son was born to Hannah Newton, young 
widow of the farmer Isaac Newton, who died a few 
months before the birth of his son. Historians and 
biographers have tried in vain to find any person of 
unusual brain power or unusually strong personality in 
the ancestry of this child. His father was ‘‘a wild, ex- 
travagant and weak man’’—at least according to testi- 
mony of his stepfather, the Reverend Barnabas Smith, 
who married Hannah Newton a few years after the 
death of her first husband. The boy himself was pre- 
maturely born and so weak tliat he was not expected to 
survive. However, he lived and became one of the 
greatest scientific minds of all time, the man of whom 
Alexander Pope wrote: 


“Nature and Nature’s laws lay hid in night. 
God said, Let Newton be!—and all was light.” 


Mathematicians, physicists and astronomers claim 
Isaac Newton as their own. However, engineers also 
have a good claim to a certain portion of his scientific 
interests and creative ideas. Back of Queen’s College 
in Cambridge, England, an old wooden bridge spans the 
Cam River. The story goes on the Cambridge campus 
that the stress analysis of this bridge was worked out by 
Newton, who was Lucasian Professor of Mathematics 
at the University. The influence of Newton’s concep- 
tions of dynamics was so far reaching that it would be 
difficult to follow up all of the contributions of the New- 
tonian Mechanics to engineering and engineering sci- 
ence. However, aeronautical engineers might be inter- 
ested in his contributions to fluid mechanics and 
especially to the problem of air resistance. 


LAws OF FLu1D MorION 


The first systematic treatment of the equations of 
fluid motion was given by Euler about a half century 
after Newton’s work. Newton’s methods were specu- 
lative and empirical rather than strictly logical and 
analytical. For example, the first edition of his 
Mathematical Principles of Natural Philosophy contains 
the erroneous statement that the velocity of a fluid 
issuing through a small orifice in the bottom of a vessel 
filled to the height 4 is equal to ~/ gh, although Toricelli 
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before him recognized the fact that the velocity of efflux 
through an orifice is equal to ~/ 2gh, i.e., the velocity of 
free fall of a body starting with zero velocity at the 
height of the water level. Newton was led to his 
erroneous statement by observation of the fluid volume 
leaving the orifice in unit time, which due to the con- 
traction is considerably smaller than the product of the 
cross-sectional area of the orifice and the velocity +/ 2gh. 
One of his collaborators called his attention to the error, 
and, reconsidering the case, Newton discovered the vena 
contracta and determined by an approximative theo- 
retic calculation the contraction coefficient as +/ 1/2, 
which is not too far from the actual value. : 


AIR RESISTANCE 


Newton’s formula for the resistance encountered by a 
body moving uniformly through a fluid medium was the 
standard formula for computation of air resistance and 
wind pressure on structures through more than two 
centuries. The development of our knowledge of air 
resistance is certainly one of the interesting chapters in 
the history of science. 

The ancient philosophers—the term philosopher in- 
cluded at that time any scholar of natural sciences—did 
not have the true concept of air resistance. Asa matter 
of fact, the Aristotelian philosophy and science had the 
vague idea that the air exerted some force on the pro- 
jectile which kept it moving. Aristotle himself de- 
scribed this motion of a body in air in a rather interest- 
ing way. He believed there could not be any vacuum 
in nature, so that the vacuum that was created by the 
displacement of the body was filled by the air streaming 
into the evacuated space and the inflowing air pushed 
the body forward. Leonardo da Vinci, who was less of 
a philosopher than an artist and a keen observer of 
nature, recognized the existence of fluid resistance and 
also the vortex formation behind a body as the main 
source of its resistance. As a matter of fact, he in- 
vented streamlining as a method of reduction of resist- 
ance long before science justified it. 

Galilei, who, for some unknown reason, is called by 
his first name, Galileo, was the first scientist who clearly 
separated idealized motion without friction and resist- 
ance from actual motion and clearly stated that the 
effect of the medium is not that of pushing a body but 
decelerating it if it moves without action of other forces, 
or decreasing its acceleration if it moves under action of 
gravity. He made experiments to find the relation be- 
tween air resistance and velocity; he found in the case 
of pendulums that ‘‘they have been resisted by air pro- 
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portionally to their velocities,” which is true if the 
resistance is produced merely by the viscosity of the 
medium. However, Galileo did not attempt to deduce 
the law of air resistance from the laws of fluid motion. 


NEWTON’S DEDUCTION 


Newton’s deduction of the law of air resistance is 
based on correct mechanical principles under the as- 
sumptions made by him. He assumed that the major 
part of fluid resistance is due to the inertia of the fluid 
matter. Under this assumption he found that the 
resistance is proportional to the density of the fluid, to 
the square of the velocity of the body relative to the 
fluid and the square of the linear dimensions of the body. 
He clearly stated the law of mechanical similarity 
corresponding to his assumptions; therefore, he has to 
be credited with the discovery of the drag formula which 
is used every day in aerodynamic calculations. He then 
went further and tried to compute the magnitude of the 
fluid resistance for bodies of a given shape and also the 
optimum ogive, i.e., the shape for the minimum resist- 
ance of a projectile. In these calculations he assumed 
that the fluid consists of equal and equidistant particles 
that have no mutual actions whatever. This assump- 
tion led him to his erroneous conclusions: first, that 
the air resistance is determined merely by the shape of 
the frontal part of the body, because only this part is hit 
directly by the particles of the fluid; second, that the 
contribution of an arbitrary surface element depends on 
its inclination only, being proportional to the square of 
the sine of the angle of inclination between the element 
and the direction of motion. 

It has often been stated, and it is true to a certain 
extent, that the common belief in the correctness of 
Newton’s theory of air resistance was an impediment in 
the solution of the problem of mechanical flight. In 
fact, the strict application of Newton’s theory gave a 
pessimistic prognosis for the feasibility of practical fly- 
ing machines. According to this theory the force acting 
on an inclined surface is produced by inelastic impact of 
the mass which hits the surface in unit time. This mass 
depends on the frontal area S sin a and is equal to puS 
sin a, where p is the density of the air; S, the surface 
area; v, the relative velocity between the air stream and 
the surface; and a, the angle of inclination (angle of 
attack). The change of momentum of the unit air 
mass is equal to the normal component of the velocity, 
i.e., tovsin a. Therefore, the normal force is given by 

N = pSv? sin? a (1) 
Thus the lift is equal to 
L = Ncos a = pSv" sin? a cos a 
and the drag 
D = Nsin a = pSv’ sin? a 


Hence the lift and drag coefficients would be given by 
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2 sin? a cos a 


(2) 


Cr = 
and 
(3) 


It is evident that according to this theory fair values of 
C,/Cp can be obtained only for small values of a; 
however, for such values the lift coefficient C, would be 
extremely small. For example, if it is assumed that 
C,/Cp = 10, Cp = 0.02. The effect of skin friction 
being neglected in this calculation makes the conditions 
even worse. 

The theoretic lift coefficient of a plane airfoil of in- 
finite aspect ratio, according to modern theory, is equal 
to 


Cp = 2sin' a 


C, = 2rsin a 


If it is assumed that the average vertical momentum 
transferred to a unit air mass involved in the lift pro- 
ducing process is equal to v sin a, the mass of air in- 
volved in the process, according to modern theory, is 
equal to pucr per unit span, where c is the chord. This 
mass corresponds to an air stream of the width 3.14c, 
whereas the mass involved in the same process, accord- 
ing to the Newton theory, corresponds to an air stream 
of the width ¢ sin a. This explains the tremendous 
discrepancy as far as the magnitude of the lift is con- 
cerned. 

The discrepancy was discovered by experiments soon 
after the publication of the Mathematical Principles, 
and the writer does not believe that the progress of flight 
was really delayed by Newton’s authority. However, 
Newton’s formula for the air pressure on inclined sur- 
faces was repeated in hundreds of books and official 
specifications. It is believed that the civil engineers 
were the last ones who recognized the discrepancy be- 
tween the old theory and the experimental facts. In 
building codes of several countries, and states and coun- 
ties the wind pressure on inclined roofs was specified 
according to Newton’s formula as late as the last dec- 
ade. This is really a remarkable proof of the inertia of 
official specifications, since, according to experimental 
evidence and also the modern theory, the wind might 
exert a lifting force on a roof consistitig of two slightly 
inclined surfaces, whereas Newton’s theory yields a 
downward directed force. 

Newton himself was dubious as to how far his results 
obtained for an idealized fluid medium could be applied 
to real fluids. It seems he considered the air as “a 
medium consisting of discrete particles.” Then he 
says: ‘In continued medium, as water, hot oil, and 
quicksilver, the obstacle as it passes through them does 
not immediately strike against all the particles of the 
fluid that generated the resistance made to it, but 
presses only the particles that lie next to it, which press 
the particles beyond, which press other particles, and 
so on; and in these mediums, the resistance is dimin- 
ished by one-half.” 

(Continued on page 548) 











Proportioning a “Canard” Airplane for 
Longitudinal Stability and Safety 
Against Stall 


JOSEPH V. FOA* 
Unwersity of Minnesota 


SUMMARY 


The conditions for static longitudinal stability and controllabil- 
ity and for safety against stall of a canard-type airplane are in- 


vestigated. 

It is found that the designer’s freedom of initiative in propor- 
tioning a canard is restricted in very definite ways. The fact 
that some of these restrictions have probably not been entirely 
known may serve to explain some of the failures of the past. 

The rules that must be followed in order to secure longitudinal 
stability and controllability and safety against stall are deter- 
mined for both the conditions of power-off and power-on, whether 
the high-lift devices are operating or not. 

High-lift devices on the tail only, as they have been used in two 
recent canards, are found to be dangerous, unless the range of the 
allowable c.g. locations is greatly reduced. The use of tail high- 
lift devices when the wing is equipped with flaps is recommended, 
but it is found that the wing and tail high-lift devices must be 
operated in conjunction with each other and under a single con- 
trol. The relation that must exist between their deflections is 


established. 
The procedures for the selection or determination of the airfoil 


section, aspect ratio and area of the tail surface, tail length and 
angle of setting of the stabilizer are indicated. 
It is felt that the tail-first airplane still deserves the earnest 


attention of designers. 


INTRODUCTION 


tie SEVERAL ADVANTAGES that the tail-first or 
canard-type airplane offers over the conventional 
airplane have attracted the attention of designers since 
the earliest days of aviation. However, to date, prac- 
tically all attempts toward a satisfactory design of a 
canard have failed. The consequence of this dis- 
couraging record has been the total disappearance of it 
in recent years. 

The proportioning of a canard for longitudinal stabil- 
ity and controllability and for safety against stall is 
quite different from that of a conventional airplane. 
The occurrence of a spin in a canard would be extremely 
dangerous; therefore, a canard must, and fortunately 
can, be made safe against stall. On the other hand, 
such safety cannot be achieved without a sacrifice of the 
controllability and of the range of the allowable c.g. 
locations. In order to attain all the advantages the 
canard can offer, it is necessary to observe a number of 
strict rules so that little freedom is left to the designer 
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of a canard in proportioning his plane. For instance, 
the effectiveness of the tail must be kept within two 
usually close limits: a maximum limit, which is imposed 
by the condition of static longitudinal stability, and a 
minimum limit, which is imposed by the condition of 
longitudinal controllability. Analogous restrictions are 
found in the range of the allowable c.g. locations, the 
selection of the aspect ratio and airfoil section of the tail 
surface, the operation of the high-lift devices, and so on. 
If such restrictions are not known and observed, a fail- 
ure is almost inevitable. 

It is the purpose of this paper to determine some of 
the above-mentioned limitations, thereby offering a 
probable explanation for some of the failures of the 
past and suggesting a possible remedy for future de- 


signs. 
SYMBOLS 


= aerodynamic center of the wing 
mean aerodynamic chord of the wing 
area of wing 
area of tail surface 
= total pitching moment about the c.g. (positive when 
stalling) 
= chordwise force on wing 
lift of wing 
lift of tail 
L,/ (p/ 2)S; ¥? 
pitching moment of wing about the A.C. 
dC,/da 
dC1,/da 
number of propellers 


Fig. 1 shows the meaning of other symbols as used in 
this paper. h,/ andj are measured from the A.C., posi- 
tive when forward; k is measured from the wing chord, 
positive when the c.g. is above the chord. All other 
symbols have the standard meaning. 

The normal force component on the wing is assumed 
to be equal to the wing lift. 
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The direction of the resultant air reaction on the tail 
is assumed to be perpendicular to the wing chord and 


equal to the tail lift. The location of the center of pres- 
sure on the tail is assumed to be constant. 

The pitching moments contributed by the body, 
nacelles, fins and rudders and by the propeller thrust are 


neglected. 
The chordwise force coefficient C, is approximately 
C. = aC, — Cp 
Write 
Cp = Co, + Cr, 


where Cp, is the parasite drag coefficient, which is here 
assumed to be constant, and Cp; is the induced drag 
coefficient: 


Co; => a;Cr 


where a; is the induced angle of attack. 

Let a. be the angle of attack of the wing measured 
from zero lift, which would give lift coefficient C, with 
two-dimensional flow about the wing profile, and let 
—ay be the geometric angle of attack of zero lift. 


Then 
a— & = As — AH 
whence 
C. = (ae — a)C, — Cr, 
But 
a. = C,/ma 


where m., is the slope of the lift curve of the wing with 
infinite aspect ratio. Then 


C. = (Cr?/m ao) — aC, — Cry 
and 
dC./dC, = 2a,, — a 
The down-wash effect of the tail on the wing is neg- 
lected. Therefore 
dCx,/dC, = (dCr,/da)/(dC,/da) = m,/m 


The lifting force acting on the plane of rotation of 
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each propeller when its axis is inclined to the direction 
of motion is approximately* ® 6 
F = CrpSV*a,/2 
with 
= |S ~ 2] 


ee sls ~ sar 


where D is the diameter of the propeller, Cg is the torque 
coefficient defined as Torque/pn*D*, J is the advance- 
diameter ratio V/nD, a, is the inclination of the axis to 
the direction of motion, and K is a coefficient depending 
on J (Fig. 2).* 


StaTic LONGITUDINAL STABILITY AND EQUILIBRIUM 


The total moment about the c.g. may be written in 
the form: 


M = L,(l — h)c + My + Che — Lhe + 3) F(j — h)c 
Dividing throughout by pV?Sc/2, 


S 
Cn = Coi(l — h) + (Cu) + 


Cy? 
(i — Co, k — Gh + DCr(j — hay (1) 


M on 


The condition for equilibrium at the incidence at 
which the wing lift coefficient is C, is C,, = 0, whence 


Oe” oe 


Cogk — (Caw — DCe(5 — Mas | (2) 


| cut + ak — ak) + 


The angle of setting of the tail, or the equivalent ele- 
vator deflection, for equilibrium at any given angle of 
attack of flight and for any given c.g. location can be 
obtained from Eq. (2). The angle of setting of the tail 
(angle between the chord of the tail and the chord of the 
wing) is the difference between the angle of attack of the 
tail, which gives the tail lift coefficient Cz,, and the 
angle of attack of the wing, which gives the wing lift 
coefficient C;,. 

From Eq. (1), 


Rew + te — 04 
m S 


1 , 
m Cri — Hh) ©) 


For stability, dC,,/dC, must be negative. 
Eq. (3) shows that if k is positive (c.g. above chord) 
the stability decreases as the angle of attack increases; 


* New information on this propeller force has been published 
recently by Rumph, L. B., White, R. J., and Grumman, H. R., 
Propeller Forces Due to Yaw and Their Effect on Airplane Stability, 
Journal of the Aeronautical Sciences, Vol. 9, No. 12, pp. 465-470, 
October, 1942. ; 
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if k is negative (c.g. below chord), the stability increases 
with the angle of attack. It also shows that the sta- 
bility is greater with power-on if j — h is negative (plane 
of rotation behind the c.g.) and greater with power-off 
if 7 — h is positive. 

For neutral stability, dC,,/dC, must be zero, whence 


Si 


mi 1+ (2a. — ao)mk + > Cri 


(4) 





ho = 


S 
m+ m+ Cr 


Eq. (4) gives the locus (io, k) of neutral c.g. loca- 
tions, i.e., the locus of the most rearward allowable 
c.g. locations for stability. For k > 0, the critical value 
of a, is that corresponding to the maximum lift coef- 
ficient CLmar.; for k < 0, io can be calculated with a, = 
0. 


CONTROLLABILITY AND SAFETY AGAINST STALL 


Eq. (2) shows how, for a given angle of attack of equi- 
librium, the required lift coefficient increases with the 
distance of the c.g. from the aerodynamic center of the 
wing (for flight test confirmation see reference 2). For 
every angle of attack of flight in equilibrium, there exists 
a limit forward c.g. location, which is that for which 
the required Cz,, is the maximum lift coefficient of the 
tail, Cr;,; and, with a given Cz,, the maximum angle 
of attack at which the airplane can fly in equilibrium 
decreases as the c.g. moves forward. 

It follows that the airplane is safe against stall if the 
tail lift coefficient required for equilibrium at the stall- 
ing angle of attack of the wing is Cr,, when the c.g. is in 
its limit rearward location (hk = h,). 

Conversely, the condition that the tail and the wing 
reach the stall simultaneously determines the limit 
rearward c.g. location. From Eq. (2), assuming that 
the vertical location of the c.g. be constant, 


S , 
Cx, / - + kCes + (Cn)w + LCriae, 
(5) 





h, = r 
Cruse. + Ct, . + Cray, 


where C:, is the chordwise force coefficient of the wing 
at stall and ap, is the inclination of the propeller axis 
when the airplane is at the stalling incidence. 

If the vertical displacements of the c.g. are not negli- 
gible, Eq. (5) gives the locus (h,, k) of the most rear- 
ward allowable c.g. locations for safety against stall. 

Since the maximum incidence at which the airplane 
can fly decreases when the c.g. moves forward and since 
there is a lowest maximum lift coefficient Cz, which 
the wing must be able to reach regardless of the c.g. 
location (e.g., the lift coefficient corresponding to the 
angle of attack of best rate of climb), it follows that the 
most forward allowable c.g. location, /,, ts determined 
by the condition that the airplane be able to fly at 
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the incidence corresponding to Czy when h = hy. 
From Eq. (2), 


S 
Cx, = aa kCcs + (Cr) w + LeCrjany 
(6) 





hy = 
S; 
CLy + CL, Ss + Crap, 


where Cz, is the chordwise force coefficient of the wing 
when C,; = Cz, and ap, is the corresponding inclination 
of the propeller axis. 

If the vertical displacements of the c.g. are not negli- 
gible, Eq. (6) gives the locus (h,, k) of the most forward 
allowable c.g. locations. 


POWER-OFF AND POWER-ON 
The most rearward allowable c.g. location for static 


longitudinal stability, power-off, is: 


S 
mM, =! + (2a, 


— ao)mk 





(ho)o = (4’) 


t 
m+ m, 3 

If j is negative, or positive but smaller than (/y)o (i.e., 
if the plane of the propeller is behind the most rear- 
ward allowable c.g. location for stability with power- 
off), then the introduction of the propeller terms as in 
Eq. (4) would yield io < (ho)o. In this case, the critical 
condition is obviously that of power-off, and (Ho) must 
be taken as the most rearward allowable c.g. location 
for stability. 

If j is positive and greater than (ho)o, the critical 
condition is that of power-on, and the propeller terms 
have to be considered in the computation of the most 
rearward allowable c.g. location, as in Eq. (4). 

If there are two or more propellers whose planes of 
rotation are at different distances from the aerodynamic 
center, then for the calculation of io only the propeller 
terms for those propellers whose j is positive and greater 
than (/o)o will have to be considered in Eq. (4). 

The following is found by analogous considerations. 

The most rearward allowable c.g. location for safety 
against stall, power-off, is: 


Cx, 1 ms + kCe, + (Cr) w 
3 S , 
(5’) 





(hr)o = 5 
CL, 5 + Cruse. 


In the calculation of h,, only the propeller terms per- 
taining to those propellers whose j is positive and greater 
than (h,)o will have to be considered in Eq. (5). 

The most forward allowable c.g. location for controll- 
ability, power-off, is: 


Cry h 2 + Coy + (Cu) 
‘8 S . 
(6 





(Ayo = 
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If j is negative, or positive but smaller than (/,)o, then 
the introduction of the propeller terms as in Eq. (6) 
would yield h,< (h,)o; therefore, the critical condition 
in this case is that of power-on. 

If j is positive and greater than (h,)o, the condition of 
power-off is critical, and (/,)o is the most forward allow- 
able c.g. location. 

If there are two or more propellers whose planes of 
rotation are at different distances from the aerodynamic 
center, only the propeller terms for those propellers 
whose j is negative, or positive and smaller than (/y)o, 
have to be taken into consideration in Eq. (6) for the 
determination of the most forward allowable c.g. loca- 
tion. 

Of all possible conditions of power-on, the most criti- 
cal is that of full throttle at the minimum speed. In 
such a condition J is a minimum; a,, K and Cg are 
great; and, therefore, F is at its maximum. 


Tue Use or HicH-Lirt DEVICES 


From Eq. (2), it appears that if Cr,;, were increased 


the airplane could be controlled to higher angles of 
attack when the c.g. is in a forward location. 

Tail high-lift devices have been used for that purpose 
in both the Focke-Wulf and the Stefanutti canards. 
However, it is clear that the improved controllability is 
acquired at the expense of the safety against stall unless 
the most rearward allowable c.g. location is moved for- 
ward or unless a safety device is used to prevent the 
operation of the tail high-lift devices when the c.g. is 
aft of a certain location. Such a solution is evidently 
impractical.* 

On the other hand, a tail high-lift device may be 
necessary when the wing is equipped with high-lift 
devices. The operation of the latter usually introduces 
a negative increment of (C,)w, whereby both h, and h, 
are reduced. Therefore, the condition of flaps neutral 
is critical for h,, while the condition of flaps in operation 
is critical for hy. 
tendency to reduce the range of the allowable c.g. loca- 
tions. Such a reduction can be avoided by also equip- 
ping the tail with a high-lift device. But in such a case 
the high-lift devices of wing and tail must be intercon- 
nected: Eqs. (5) and (6) give the relation that must 
exist between the maximum lift coefficients of tail and 
wing, and therefore between the respective flap deflec- 
tions, in order to keep hy and h, unchanged. Of course, 
Cr, will change when the high-lift devices are in opera- 
tion. 

With Eq. (4), 4o must be calculated for different con- 
ditions of operation of the high-lift devices, including 
the neutral condition. The greatest obtained value 
must be assumed as the most rearward allowable c.g. 


location for stability. 
* The Focke-Wulf FW 19a did not have any safety device, and 


its range of allowable c.g. locations was reduced to 3.65 per cent 
of the mean chord—i.e., only about 3.6 in. (see reference 2). 
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FURTHER LIMITATIONS 


If h, < ho, then for h, < h < ho the airplane would still 
be safe against stall but would be longitudinally un- 
stable. Therefore, in such a case 4) must be taken as 
the limit rearward c.g. location. However, that 
means a cut in the range of the c.g. locations. In addi- 
tion, as the c.g. cannot be allowed to reach the position 
h,, the airplane can never fly at the angle of attack of 
maximum lift. Hence, the condition h, < ho is undesir- 
able. 

On the other hand, if h, > io, the airplane may be too 
stable. As in the case of the conventional airplane, a 
condition near neutral stability when the c.g. is in the 


most rearward position is the most desirable. Such a 
condition is obviously achieved when 
h, = ho (7) 


If a certain range rc of c.g. locations is desired, also the 
equation 


(8) 


h;-h,=r 


must be satisfied. 

It is obvious that the conditions in Eqs. (7) and (8) 
are further limitations to the designer’s freedom in pro- 
portioning the canard. 


SAMPLE CALCULATION 


Data: 
Gross weight............. 4,000 Ibs. 
| BOR ECE OOO 200 sq.ft. 
Aspect ratio of wing...... 10 
Aspect ratio of tail....... 4 
atkey sink ik Mitechio aA atin 0.085 
ate ee ene rT 0.070 
| OPES AP eT ae rer ee 5 ft. 
eT Se dre Aes lag Aes cis 1.6 
A dicho sarma sees eb i 1.3 
Ri lt deities nih tia awake 0 
iin Sir donenndisee tik 0 10 deg. = 0.175 rad 
EN ae ae 0 
En 1 
esheets Gi tats wings Sia —1 
ge: a eee it ane ee 8 ft. 
J in the critical condition of 
ee rr 0.6 

E m ed Paphuaaela 0.025 

J 2dJ 
Ge We Fs ws we 4.8 

Padeeeevwekesaasdivasen 0.54 
Pear esha ds sees oe eee 6 in. 


Determine tail length, tail area, aerodynamic charac- 
teristics of tail profile, angle of setting of the tail for dif- 
ferent speeds and c.g. locations, minimum speed for 
different c.g. locations. From Eqs. (7) and (8), 4, = 
ho = 0.54 and hy = h, +r = 0.64. 
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Since j is negative, 4p and h, must be calculated for 
the condition of power-off. From Eq. (5’), 


S, 
CL es 300 


+ 1.6 





h, = 0.54 = 
os = 
The condition of power-on is critical for hy. 
Eq. (6), 
h, = 0.64 = 


Ss 
Ogee 4.8 X 5 X 0.025 X 0.175 
(b) 





64 
13+ Cy, Ot + 4.8 X SE x 0.025 X 0.175 


Eq. (7) yields 
0.070 
0.085 
0.070 S, 
0.085 200 

whence Cz,, = 1.32—i.e., the tail must be such as to 

reach the stall when its lift coefficient is 1.32. 

Then, from Eqs. (a) and (b), S; = 35 sq.ft., 7 = 4.27, 
and lc = 4.27 X 5 = 21.4 ft. 

The tail lift coefficient Cz, for equilibrium at different 
angles of attack and for different c.g. locations can be 
calculated from Eq. (2). The angle of setting of the tail 
is the difference between the angle of attack of the tail, 
which gives the tail lift coefficient Cz,, and the angle of 
attack of the wing, which gives the corresponding wing 
lift coefficient C;. 

In Fig. 3 the values of Cx, for equilibrium are plotted 
vs. the corresponding air speeds at sea level, instead of 
vs. the corresponding values of C,. The weight of the 
airplane is assumed to be the same for all c.g. locations. 

Fig. 3 also shows how the minimum speed increases 
as the c.g. moves forward. 


Ct,1(S;/200) 


1(.S,/200) 





CONCLUSIONS 


(1) The most rearward allowable c.g. location is de- 
termined by two conditions: static longitudinal stabil- 
ity and safety against stall. The condition of static 
longitudinal stability yields ho (Eq. 4); the condition 
of safety against stall yields h, (Eq. 5). The greater of 
the two values must be taken as the most rearward 
allowable c.g. location. 

(2) The maximum angle of attack at which a canard 
can fly in equilibrium decreases as the c.g. moves for- 
ward. 

(3) The condition of safety against stall is attained 
at the expense of the maximum lift coefficient and, 
hence, of the minimum speed at which the airplane can 
fly in equilibrium. If h, = ho so that h, is the most 
rearward allowable c.g. location, the airplane can reach, 
but not exceed, the angle of attack of maximum lift 
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coefficient when the c.g. is in the most rearward loca- 
tion. If ho > h, so that io is the most rearward allow- 
able c.g. location, then the canard can never reach the 
angle of attack of maximum lift coefficient, even when 
the c.g. is in the most rearward position. 

(4) The lowest maximum lift coefficient Cz, which 
the wing must be able to reach, regardless of the c.g. 
location, determines the most forward c.g. location 
(Eq. 6). 

(5) Both h, and h,, and therefore also the range of 
the allowable c.g. locations if h, 2 ho, are linearly pro- 
portional to the tail length. 

(6) The condition h, = his desirable. This condi- 
tion can be achieved by a proper selection of the profile 
or aspect ratio of the tail. 

(7) Since the maximum lift coefficient that the wing 
can reach decreases as the c.g. moves forward, it is ad- 
visable that the useful load be distributed so that the 
c.g. is in the most rearward location in the full load con- 
dition and moves forward as the load decreases. 

(8) High-lift devices on the tail, to be used only 
when the c.g. is in a forward location for the purpose of 
increasing Cr;—i.e., the controllability—are dangerous 
and cannot be recommended. 

(9) When the wing is equipped with high-lift de- 
vices, the use of high-lift devices on the tail, in conjunc- 
tion with those operating on the wing, is recommended. 
The wing high-lift devices alone would usually decrease 
hy; and h,, while the tail high-lift devices alone would 
increase both h, and h,. The operation of the wing and 
tail high-lift devices in conjunction with each other and 
under a single control can be done in such a way as to 
leave h; and h, unaffected. For every deflection of the 
wing high-lift devices, Crime, and (C,)w are known; 
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then, Eqs. (5) and (6) will yield Cz,, (i.e., the corre- 
sponding deflection of the tail high-lift device) and Cz, 
(i.e., the maximum wing lift coefficient that can be 
reached with the assumed deflection of the high-lift de- 
vices when the c.g. is in the most forward location). 

(10) If the wing and tail high-lift devices are not 
operated in conjunction with each other under a single 
control and according to the rules described above, then 
the range of the allowable c.g. locations is reduced; hy 
and h, have to be calculated (Eqs. (5) and (6)) for all pos- 
sible combinations of operation of the wing and tail 
high-lift devices; and the smallest of all values of h, 
and the greatest of all values of h, have to be taken as 
the most forward and the most rearward allowable c.g. 
locations, respectively. 

(11) If the planes of rotation of the propellers are 
behind the aerodynamic center of the wing (j < 0), the 
condition of power-on is critical for h, and that of power- 
off is critical for h,. If the planes of rotation of all or 
some of the propellers are located before the aerody- 
namic center (j > 0), 4, and h, must first be calculated 
for the condition of power-off [(/,)o and (h,)o, Eqs. (5’) 
and (6’)]. Then, in the calculation of h, (Eq. 6), only 
the propeller terms pertaining to those propellers whose 
jis less than the calculated (h,)o must be considered; in 
the calculation of h, (Eq. 5), only the propeller terms 
pertaining to those propellers whose 7 is greater than the 
calculated (h,)o have to be considered. ho must be 
checked for the conditions of power-off and power-on in 
the same way as h,. 
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(12) The angle of setting of the tail surface or the 
equivalent elevator deflection for equilibrium at a given 
angle of attack at which the wing lift coefficient is C, 
can be calculated by means of Eq. (2) for any c.g. loca- 
tion. Eq. (2) yields Cz, and therefore the angle of 
attack of the tail or the elevator deflection necessary to 
keep the airplane in equilibrium at the chosen angle of 
attack. 
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Application of Maclaurin Series to Beams 
Under Simultaneous Transverse and 


Axial Loading 


A. G. STRANDHAGEN* 
Carnegie Institute of Technology 


INTRODUCTION 


The method of determining the elastic line of columns and 
struts is, to start with, the fundamental equation (d*y/dx?) = 
—M/EI, and by two quadratures the deflection curve y = f(x) 
is obtained. Most of the current methods carry this procedure 
through in some manner or other. While these methods are 
practical, it can be shown that in many cases other methods are 
shorter. As an illustration, application of the Maclaurin 
formula in solving ior the deflection curve of beams was shown 
by Hetényi.! In this paper the writer extends the Maclaurin 
formula, which developed a function of real variable, say x, 
into a polynomial plus a remainder, to an infinite series known 
as the Maclaurin series, which expresses the function as a power 
series provided the remainder after m terms converges toward 
0 as n become infinite. 

Use of the Maclaurin series simplifies the solution of many 
problems and therefore should be of interest to engineers. 

There is no approximation involved in this procedure since all 
terms are included, as will be evident from the expansion in the 
illustrative examples. 





N PROBLEMS RELATING to the bending of beams under 

simultaneous transverse and axial loading, this 

method begins with the relationship between curvature 
and moment in the familiar equation 


(d2y/dx?) = —M/EI (1) 


where x and y are Cartesian coordinates, M equals the 
bending moment at the cross section through the point 
in question, J equals the moment of inertia of the cross 
section with respect to the neutral axis, and E equals 
the modulus of elasticity. The usual assumptions of 
small deflections, plane of cross section before bending 
remaining plane after bending, and homogeneity are 
implied. 

If E and J are constant, and M, the bending moment, 
varies as some function of x and y, then Eq. (1) can be 
expressed symbolically as 


F[(d*y/dx*), y, x] = 0 (2) 


a nonhomogeneous, second order, linear differential 
equation with constant coefficients, the solution of 
which can be found. 


Received May 16, 1942. 
* Instructor in Mechanics. 
t By standard procedures is meant the more common methods 
such as reduction by differentiation to homogeneous form 
operator method, variation of parameters, etc. 


y=fO+Ox+ YD 


» Spee 


Instead of solving Eq. (2) by the standard pro- 
cedures{ found in texts on differential equations, the 
deflection curve y = f(x) can be in a sense generated 
from the one equation—i.e., f’(x) = (d*y/dx*) = 
—M/EI. 

It is recalled from the calculus that if f(x) and all its 
derivatives exist and are continuous over an interval 
about x = 0 then f(x) can be expanded in that interval 
in an infinite series of the form 


y = f(x) = £0) + £’(O)x + ("(0)x*/2!) + 
(£’’"(0)x8/3!) + ... + (£"(O)x"/n!l) +... (3) 


where n increases indefinitely and f’(0), £”(0), f’’’(0), 
... £"(0) are, respectively, the first, second, third, and 
nth derivative of f(x) when x = 0. Maclaurin’s 
infinite series, Eq. (3), represents f(x) provided 
(f"(0)x"/n!) — 0 as n — o. The Maclaurin series 
expansion of f(x) is known to be unique.{ Thus 
given f"(x) = (d*y/dx*) = —M/EI, then by successive 
differentiation of the same with respect to x there 
results f’’’(x), f'Y(«), ... f(x). Evaluating f’(x), 
f’’’(x), ...f"(m) at x = O and then substituting in 
Eq. (3), a power series in x results, which in the cases 
discussed can be converted over into the circular and 
hyperbolic functions. f(0) and f’(0) are determined 
from the boundary conditions. 

Eq. (3) can be written more compactly and this form 
will be found convenient although not necessary.** 


[£"(0)x"/m!] + 


y [f"(O)x"/m!] (4) 
Closely allied to the Maclaurin infinite series is Taylor’s 
infinite series, which is used in some of the solutions 
of problems and is as follows: 


y=(O+*Oe-a2+ FE x 


n=2,4,... 


[f"(a) (x — a)"/n!] + 2 [f"(a)(x — a)"/n!] (5) 


which states that the function y = f(x) is evaluated in 


} Derivation of Maclaurin and Taylor series can be found in 
Wilson’s Advanced Calculus. 

** The rearrangements of the terms of the power series are 
valid within the interval of convergence. See, for example, 
Wilson’s Advanced Calculus, p. 442. 
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terms of the value of function and its derivatives at 
x = ain powers of (x — a). 

The following expansion in the power series of func- 
tions encountered frequently in this paper are given 
here without proof.* 


nu=@ 


p> 


n=2,4,... 





(—1)"2k"x"/n! (6) 


(cos kx — 1) 


n= @ 


» 


n=3,5,... 


(sin kx — kx) (—1)%-0/2p%"/y! (7) 


ae @ 


> 


n=2,4,... 


(cosh kx — 1) k"x"/n! (8) 


‘Sy - 


2=3,5,... 


(sinh kx — kx) k"x"/n! (9) 


APPLICATIONS 


In what follows, unless the contrary 1s specified, it 
is assumed that the moment of inertia J is constant, 
that the material is homogeneous, and that E is constant 
everywhere. It is assumed that the bending is uni- 
planar—i.e., that the bent centerline of the column 
lies entirely in one plane. The condition for this is 
that the plane of bending contain one of the principal 
axes of the cross section. The load P is applied at 
the center of area of the cross section and along the 
unstrained centerline. For convenience, let k? 
P/EI in the following illustrations. 


Case 1. Cantilever Beam Under Axial Compression and 
Transverse End Load (Fig. 1) 


At any point x 
(d*y/dx*) = [QU — x)/EI] + [P(@—y)/EI] (1.1) 


From the end conditions: The deflection is zero when 
x = 0, the slope to the deflection curve is zero when 
x = 0, and when x l the deflection is equal to 6. 
Thus, the boundary conditions may be expressed 
as 


{0)=0 f'(0)=0 fd) =6 


Then by successive differentiation of Eq. (1.1) with 
respect to x and expansion about x 0 and by the 
substitution of k? = P/EI, there results 


(1.2) 








i"(x) = we £"(0) = +R 
ee) = —S rm) = - 2 
Me) = gM) = — Se es 
Mm =-2- M0) =+2e 


The nth term of f"(0) is 
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Fic. 1. 
(0) = —(—1)”2(Ql/EDk"-? — (—1)™*k"6 
where m = 2, 4,6, ... 
{"(0) = (—1)-/2(0/EDk"-* where n = 3, 5, 7, ... 
From Eq. (1.2), f(0) = f’(0) = 0 and employing f*(0) 
from above, substitution in the Maclaurin infinite 


series, Eq. (4), gives 


n= @ 


pees 


"S$ (—1e%x"/nl] + Q/ED > (—1)0-v/a 
n=2,4,... n= Rica 
(k"—8x"/n!) 


From (6) and (7) one can recognize the expansion of 
(cos kx — 1) and (sin kx — kx) after multiplying the 
first summation by k?/k? and the third by k*/k* and 
remembering that k? = P/EI, thus 


y = — (QIl/P)(cos kx — 1) — 8(cos kx — 1) + 
(Q/Pk)(sin kx — kx) (1.3) 


In Eq. (1.3) the value of 6 may be obtained by the use 


of the third boundary condition y = 6 when x = /. 
Then 

6 = (Q/P)[(tan kl/k) — 1] (1.4) 
atx =/, The maximum bending moment is /,,,,. = 


Ql + Pé = (Q/k) tan ki at x = O and the equation 
of the elastic line is obtained by substituting the value 
of 6 from Eq. (1.4) into Eq. (1.3) 


y = —(Q/Pk) tan ki(cos kx — 1) + (Q/PkR) sin 
kx — (Qx/P) 


Case 2 


Case 2 is the same problem as in Case 1, except that 
the function is to be expanded about x = /. From 
Fig. 1 the slope is f’(/) and is not zero, and f(/) = 6. 
Then 


f"() = 0 

f’’"(1) = —k*f’() — (Q/ED 
Yl) = 0 

() = 


7 + (Q/ENk? 
The nth term of f*(/) is 
i"(1) = 0 


f"(I) = (— 1) —0/26/(]) R21 + 
(- 1)%"—0/2(QR"-38/ET) when n 


when = 2, 4, 6, ... 


3,5, 7, ..- 


























APPLICATION OF MACLAURIN SERIES TO BEAMS 


In this case x = / has been arbitrarily chosen as a point 
of expansion; thus Taylor’s infinite series (Eq. (5)) is 
used where a is replaced by /. Then 


y=6+fO@-)+IOX 


7 it. 1)" —D/2p—-1(% las 1)"/n!] > 
ae 

Q — — 1(n—1)/2hn-3(,. — ]\n 

Er, a2, | 1 k™-*(x — 1)"/nl] (2.1) 


By inspection of Eq. (2.1), the odd n’s suggest a sin kx 
expansion. Multiply the first summation by k/k and 
the second by k*/k*. From Eq. (7) and factoring, 
there results 


y = 6 + (f’(2)/k) sin R(x — 1) + (Q/Pk) X 
[sin k(x — 1) — k(x —1)] (2.2) 
when x = 0, y = 0; then 
5 = —(Q/Pk)(kl — sin kl) + [£’(1)/k] sin ki (2.3) 


when x = 0, f’(0) dy/dx = 0; differentiating Eq. 
(2.2) and placing x = 0, 
f’(l) = (Q/P)[(1 — cos kl)/cos ki] (2.4) 


Substituting f’(/) from Eq. (2.4) into Eq. (2.3), 6 
becomes 


I| 


6 = (Q/P)[ (tan ki/k) — 1] 


This result is the same as Case 1, Eq. (1.4), thus demon- 
strating that the expansion may be taken about either 
end of the beam. 


(2.5) 


Case3. IdealColumn. Position Fixed at Each End (Fig. 2) 


Both ends are free in direction but fixed in position. 
This case, like many cases studied in elementary 
strength of materials, will demonstrate the ease and 
simplicity of the method. 

Take the origin at 0, the midpoint of the length of 
the beam, then 

d*y/dx? = —Py/EI 
The boundary conditions are 


f{(0) =6 f'(0) =0 f(/2)=0 








Thus 
Expansion About (0, 4) 
f"(x) = —Py/EI f"(0) = —k% 
f'"(x) = —Py’/EI f’""(0) = 0 
f'%(x) = —Py"’/EI f'Y(0) = +k48 
fY(x) = —Py""/EI fY(0) = 0 
| 
' | 
(0) = 0 when s = 3, 5,7, ... 
f"(0) = (—1)"/2k% when n = 2, 4, 6,... 
Then substituting in the Maclaurin series gives 
y=6+6 YO (-1)%2%R"x"/n! (3.1) 
n=2, 3,4... " 
or 
y = 6+ 6 (cos kx — 1) = bcos kx (3.2) 
when x = //2, y = 0, Eq. (3.2) becomes 
cos (kl/2) = Owhen ki/2 = nx/2 where n = 1,3,5,... 
ork = x/l 
k? = ?/[? = P/EI 
or P,, = r°EI/l? 


Case 4. Bending of a Compressed Bar by a Couple (Fig. 3) 


Both ends are fixed in position but free in direction. 
A couple Mz is applied at B. Take the origin at 0. 
Then 


d*y/dx* =[—(M,/l)x — Py}(1/ED) 


Boundary conditions are 





f0)=0 f) =0 (4.1) 
f"(x) = | - as = Py|t f"(0) =0 
f’"(0) = — ae — k*f’(0) 
f(x) = —Py"/EI f1%(0) = 0 
ea = 5 2 {%(0) = re +k") 
| | 
f"(0) = 0 when m = 2,4, 6, ... 


f"(0) = (—1)-0/2(M,/PI)k"-" +. 
(—1)™-/2'(0)k"—! when n = 3, 5,7, ... 


Substituting in the Maclaurin series gives 


rOneee 





» 


(- 1) @~D/2(px—14" /n!) (4.2) 


+ £'(0) 
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or 
y = f’(0)(sin kx/k) + (Mg/PIk)(sin kx — kx) (4.3) 
when x = /, y = 0; then by substitution in Eq. (4.3) 
and solving for f’(0) 
f’(0) = (Mz/P)[(R/sin kl) — (1/))] 
Substituting f’(0) from Eq. (4.4) into Eq. (4.3), then 
the equation of the elastic line is 


y = (M,/P)[(sin kx/sin kl) — (x/l)] 


(4.4) 


* (4.5) 


Case 5. Beam on End Supports Under Axial Compression 
and Triangularly Distributed Transverse Load (Fig. 4) 
At any point x 

d*y/dx? = f{"(x) = [—Py — (wlx/6) + (wx'/6l)]1/EI 

The boundary conditions are 


{(0) = 0 


(6.1) 


f(l) = 0 
Expansion of f”(x) about x = 0 can be shown to give 
the following f”(0). 
{"(0) = 0 


(0) = (—1)—-/2~"—1/(0) +. 
(—1)%-P/22"—8w)/6EI when n 


when n = 2, 4, 6,... 


£"(0) = (—1)%- /2k"-5 w/IEI when n 





Fic. 4. 


Thus, substitution in the Maclaurin series gives 


(- 1) &»—D/2pn—1yn 


y=fOx+f0) bX 
n=3,5,... n! 
wl = (—1)(*-D/2pn—an 
6EI n=3,5,... n! 
w n= (—1)(*—D/2p2—5yyn 


LEI n=5,7,... n! 











and 

y = (f’(0)sin kx/k) + (wl/6PR)(sin kx — kx) + 
(w/LEIk®)(sin kx — kx + (k*x*/6)] (6.2) 

when x = /, y = 0; then from Eq. (6.2) 


(0) = w/P[(1/k sin kl) — (1/lk*) — (l/6)] (6.3) 


Substituting f’(0) from Eq. (6.3) in Eq. (6.2) the 


equation of the elastic line becomes 


y = w/P[ (sin kx/k? sin kl) + (x?/61) — 
(/x/6) — (x/k*l)] 
REFERENCE 


1 Hetényi, Miklés, Deflection of Beams of Varying Cross Section, 
Journal of Applied Mechanics, Vol. 4, pp. A-49-A-52, June, 1937. 


















Elliptic Tunnel-Wall Corrections 
on Drag and Stall 


J. G. GAVIN, JR.,* ano R. W. HENSELt 


SUMMARY 


The tunnel-wall effects on drag and stall, given by the usual 
simple corrections, are not adequate when the span-diameter 
ratio is high or when the spanwise lift distribution varies con- 
siderably from elliptical. In order to determine more exact inter- 
ference corrections, especially for elliptic throat tunnels of closed 
section, the upwash induced by the tunnel walls on the major 
axis, due first to the presence of lifting lines of different span 
lengths and then to the presence of typical wing distributions 
built up from these lifting lines to simulate actual wings, is com- 
puted for a number of wing distributions and spans. All the cal- 
culations are for wings whose aspect ratio is eight, and all are 
based on the constants of the Wright Brothers elliptic throat wind 
tunnel at the Massachusetts Institute of Technology. 

The magnitude of these corrections is found to vary inversely 
with the taper ratio of the wing, the corrections being less for a 
highly tapered wing and more for a rectangular wing of the same 
span and aspect ratio. For those cases where the wing tips lie 
outside the foci of the throat section, the variation of the induced 
angle of attack along the span is no longer negligible at high lift 
coefficients. Its geometric effect is to cause an apparent washin 
angle toward the wing tip. The induced drag corrections have 
minimum values when the wing tips are at or near the foci of the 
tunnel for all usual lift distributions along the span. As the 
wing tips approach the tunnel wall, both the induced drag correc- 
tions and the angles of apparent washin or of apparent twist 
increase rapidly. 


SYMBOLS 

A = semimajor axis of tunnel cross section 
B = semiminor axis of tunnel cross section 
b = model wing span 
c = focal distance of tunnel cross section 
Cr = lift coefficient 
ACp; = induced drag coefficient correction due to tunnel walls 

= lift (local force along span) 
A.R. = aspect ratio 
Sr = tunnel cross-sectional area 
S» = model wing area 


chord at wing root 

fictitious tip chord 

U = undisturbed stream velocity 

induced upwash due to tunnel walls 

distance from center of tunnel cross section along A 





taper ratio = 


— 
= 


bf = 

gE = hyperbolic component in hyperbolic, ellipsoidal co- 
ordinate system 

3 = £, corresponding to a given wing tip position (tip in- 
side c) 

n = ellipsoidal component in hyperbolic, ellipsoidal co- 
ordinate system 

No = constant 7 giving ellipse representing tunnel wall 


Presented at the Aerodynamics Session, Tenth Annual Meeting, 
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a = n/no for n between 0 and mo 

a’ = a, corresponding to a given wing tip position (tip be- 
yond c) 

8 = y/(b/2) 

a = angle of attack 


Aa = induced angle of attack due to tunnel walls 

r = local circulation around wing 

To = circulation at wing root 

K = T segment/T 

50 = upwash correction factor 

im = mean upwash correction factor for a given loading and 


span 


TUNNEL CONSTANTS 
B/A = 0.732 c = 0.68A 
k = 1.584 Sr = kre? 


THe CASE OF THE LIFTING LINE IN THE ELLIPTIC 
TUNNEL 


hewn FIRST STEP is to determine the tunnel-wall inter- 
ference on a lifting line that lies symmetrically on 
the semimajor axis of the tunnel. Following cus- 
tomary procedure, the wall of the tunnel is replaced by 
a system of reflected vortices. Because the position of 
the vortices is not immediately apparent in the elliptic- 
hyperbolic coordinate system determined by the tunnel 
section, it is necessary to convert the system into a 
rectilinear system by means of the following trans- 
formation :! 


(1) 


In this rectilinear system it is possible to fix the posi- 


Z = cos (yn + #€) 






















Fic. 1. Reflected vortex pattern elliptical wind tunnel. 








-—Focal points. I—Wing tips inside foci. F—Wing tips at 


I.Ae.S., New York, January 30, 1942. 
foci. O—Wing tips outsiue foci. 


*T Honors Group Students in Aeronautical Engineering. 






533 





534 





tion of the reflected vortices so that these, together 
with the vortices of the lifting line within the tunnel, 
produce the desired boundary conditions at the tunnel 
walls. These vortex positions can then be transferred 
back into the elliptic-hyperbolic system. This result 
is shown in Fig. 1 for three cases: tips of lifting line 
inside foci; at foci; and outside foci. This vortex pat- 
tern extends to infinity. 

Using the vortex pattern a summation can be made 
to determine the wall-induced upwash at any point 
along the semimajor axis of the tunnel, either inside or 
outside the tips of the lifting line. The upwash, in the 
case where the tips of the lifting line are outside the 
foci of the tunnel, is: 


n=o@ 


yin x 


n=2 


T 
hag 





1 
ia (n — a)no — Bcoshany 


1 
cosh (% — a) + 8 cosh ano 





1 
cosh (7 + a) — 





8B cosh ano 





1 
7 
cosh (n + a)no + 8 cosh “4 @) 


where 7 is an even integer. 
Similarly, for tips of the lifting line inside the tunnel 
foci the expression is: 
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Fic. 2. Spanwise variation of upwash correction factor for 
lifting lines of varying lengths. 
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we cos I’ cos & 
orc bs at "x 


n=2 





a cosh ny + 8B + 
(cosh nyo + 8)? cos? — + sinh? mm sin? £ 


cosh mm — B (3) 
(cosh m9 — 8)* cos? £ + sinh? mm sin? & 





where is an even integer. 

These series converge rapidly and the first three terms 
yield accuracy within approximately 1 per cent. For 
convenience, these expressions can be rearranged, in- 
troducing a tunnel constant, k = 1.584. 


Tips outside foci: 


Tc cosh ano 
W = ——— 
2kac* 
Saal +1 same bracketed 
— 2 cosh ano he ae »! Boose as Eq. ae 
(4) 
Tips inside foci: 
Tc cos & 
W = 
2krc* 


— +1 ' 
b _1)? same bracketed 2 
{ tH iia Bee as Eq. (3) (5) 


where c cosh am and ¢ cos — represent the semispans of 
the lifting lines for tip positions determined by a and &. 
do is the term included within the braces in Eqs. (4) and 
(5). The variation of 5) has been calculated for tip posi- 
tions corresponding to: 


& = 0.457, 0.47, 0.32, 0.157, 0 
= 0.5, 0.75, 1.00 


For each lifting line tip position, 5) has been calculated 
for various values of 8 both on and beyond the span 
of the lifting line. The results are shown in Fig. 2. 








Lirt DISTRIBUTIONS APPROXIMATED BY ASSEMBLIES OF 
LiFTInG LINES 


The next step is to make use of these results to find 
the upwash variation along wings having normal lift 
distributions. The four lift distributions used* * are 
shown in Fig. 3a, which are distributions in the presence 
of the tunnel walls. Any given lift distribution can be 
approximated by a set of lifting lines having the correct 
length and strength. It must be remembered that it is 
necessary to take account out to the wing tip of the 
wall-induced upwash beyond the tips of all lifting 
lines that have spans shorter than that of the wing 
under consideration. In order to handle this operation, 
Eqs. (4) and (5) are modified further. From Eq. (5), 


= (Ic cos £/k2mc?) do 


(6) 
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Fic. 3a. Lift distributions used. 
Fic. 3b. Breakdown example. Elliptical distribution, tips at 
foci case. 


This may be applied to a given lifting lines segment m 
of the total lift distribution to give the upwash at a 
given point p on the span of the wing: 


G).- 


A summation of the several segments yields 
WwW T> | > | 
<a Kn€ COS £m 50m 8 
(F ) Uk2rce? Lm ( En bOme) (8) 


where k,,I) = I'm. For the elliptic distribution with 
tips inside the tunnel foci, it can be shown that 


T)/U = C,Sw/me cos ¢’ (9) 
By substituting in Eq. (8), remembering that S; = 


krc?; 
W CLSw [2 : ] 
il ee an Kin€ COS Em50m 10 
”) 2mrc cos §’Sp Lom ( 008 faite) | (10) 


This is the upwash at a given point p on the span of a 
wing whose semispan is ¢ cos £’ and whose lift distri- 
bution is approximated by m lifting lines. Similarly, 
for the case where the wing tips lie outside the foci, 


—} cs K€ COSh A_7050m 11) 
(7 » 2c cosh a’mSrLlm ( nobOm) | ( 


Then, Eqs. (10) and (11) can be generalized to form a 


single expression: 
W CiSw p | 
—)} = — Kn¥m50m, 
(7) rbSr m ( . °) 


T'mC COS Em 
Uk2ac2 


(7) 





(12) 
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where } is the wing span and y,, is the semispan of the 
lifting line m. 

The term C,Sy/2bS,7 of Eq. (12) has been developed 
for the case of an elliptic lift distribution. However, 
by assuming a constant aspect ratio, the total lifts for 
all distributions for a given wing tip position are con- 
stant at a given lift coefficient. Hence, for other dis- 
tributions, 7 is replaced by r(T0¢17,/To omner)- 

Fig. 3b shows a typical breakdown of a lift distribu- 
tion into its component lifting lines and illustrates how 
the factor K,, is determined. The lengths of the lifting 
ines are chosen to coincide with those used in the 
calculations of “The Case of the Lifting Line...” so 
that Fig. 2 can be used directly to find 50,.. To deter- 
mine the induced upwash variation along a given wing, 
the summation indicated by Eq. (1%) is performed at 
several points p along the wing span. Fig. 4 shows the 
type of induced upwash variation occurring in cases 
with the wing tips inside, at and outside the tunnel 
foci. 





WING TIPS INSIDE FOCI 





WING TIPS AT FOCI - ELLIPTICAL LOADING 





WING TIPS OUTSIDE FOCI 











Fic. 4. Typical upwash distributions. 


INDUCED UPWASH VARIATION 


For each of the lift distributions used, wings having a 
constant aspect ratio of 8 were taken with their tips at 
the following five positions: 


(1) 0.3% 
(2) 0 

(3) 0.5 y’ = 1,111C 
(4) 0.75 = 1.255C 
(5) a’ = 1.0 = 1.468C 


For each tip position the spanwise variation of the up- 
wash corresponding to the four lift distributions was 


b/2A = 0.40 
b/2A = 0.68 
b/2A = 0.76 
b/2A = 0.86 
b/2A = 1.00 


0.585C 
1.000C 
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plotted. Typical plots are shown in Figs. 5 and 6. As 
an illustrative example, the spanwise variation of wall- 
induced upwash has been reduced to the spanwise 
variation of induced angle of attack. This is shown in 


Fig. 7. 











Fic. 5. Spanwise variation of upwash due to tunnel walls. 
Tip position: a@ = 0.75; 1.255c; b/2A = 0.86; A.R. = 8. 
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Fic. 6. 
Tip position: 


Spanwise variation of upwash due to tunnel - walls. 
= 8. 
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20 


w 


&« IN DEGREES 
fe} 
TUNNEL_WALL 


| 











50 C 7SC Looc 25C ISOC 


¥ 


Zoc 


Fic. 7. Example of spanwise variation of induced angle of 
attack. Elliptical distribution. Cz, = 1.5; a@ = 0.75; 1.255c; 
b/2A = 0.86; A.R. = 8. 


WALL-INDUCED DRAG AND ANGLE CORRECTIONS 


In order to determine a mean upwash correction 
factor for the purpose of calculating induced drag cor- 
rection and the mean wall-induced angle of attack, an 
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average is taken of the upwash variations plotted. The 
development is as follows: 


W S, 
Acp: = Cr (*) = by my 
av, T 


(W/U)av. 
(Sw/Sr)Cz 


Thus, for each combination of wing tip position and lift 
distribution there is one value of 5. The variation of 
5y with wing tip position is plotted for each of the four 
lift distributions. This plot is presented in Fig. 8. 
For an infinitely small span the four cases merge and 


bu = 1/s(50)¢=4/2 


(13) 


(14) 


bu = 


6, 
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Fic. 8. Mean upwash correction factor. Variation with 
wing span. A.R. = 8. Wright Brothers Wind Tunnel, 
M.LT. 


An average induced angle of attack suitable for use 
as an angle correction due to the tunnel walls may then 
be taken as 


Aay = 5u(Sw/Sr) Cr (15) 


CONCLUSIONS 


Before stating any conclusions it might be well to 
note again the restrictions that have been assumed in 
the numerical computations. The tunnel is of closed 
section. The eccentricity of the tunnel section enters 
immediately into consideration; hence, all results are 
limited to the Wright Brother Wind Tunnel or one with 
a similar throat section (B/A = 0.732). From the 
point where assemblies of lifting lines are used to ap- 
proximate lift distributions of wings, the results are 
further limited to wings of A.R. = 8. The present cal- 
culations are in a form that facilitates application to 
other aspect ratios; a change in aspect ratio involves a 
constant percentage change in the values of induced 
upwash. 

For the Wright Brothers Wind Tunnel, or any reason- 
ably similar tunnel, the following generalities may be 
formed on the basis of the preceding calculations. 

(1) At high lift coefficients, when the wing tip lies 





ELLIPTIC TUNNEL-WALL CORRECTIONS 


outside the focal point, the variation of the induced 
angle of attack along the span is no longer negligible. 
This amounts to an apparent washin of angle which 
becomes severe as the wing tip approaches the wall. 
As a result, determinations of stall characteristics in 
the tunnel, when the wing span approaches the tunnel 
major diameter, are conservative—i.e., the wing tips 
will stall at higher angles of attack in free air. 

(2) Other things being constant, tunnel-wall inter- 
ference is less for lift distributions in which the lift is 
concentrated toward the centerline. That is, for un- 
twisted wings, those with high taper ratio have tunnel- 
induced upwash of smaller magnitude than that of 
wings with low taper ratios (Figs. 5 and 6). 
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aspect ratio, other conditions being held constant. 
This fact may be seen from consideration of Eq. 
(12). 

(4) By the use of a mean upwash correction factor, 
5y, as determined in this paper, the induced drag cor- 
rection can be calculated more accurately than pre- 
viously. Using 6, in the determination of the induced 
angle of attack correction yields a result that is an 
average of the wall-induced angle of attack. This in- 
duced angle of attack correction is suitable for wind- 
tunnel work with the exception of stall pattern tests 
and is more accurate than that used formerly. 

(5) For wings of normal lift distributions, the mean 
upwash correction factor, 5y, is a minimum when the 


(3) Tunnel-wall interference is less for wings of high wing tips are approximately at the tunnel foci (Fig. 8). 
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SUMMARY 


An analysis is made of the secondary stresses that arise in flat- 
sheet-covered box beams in which the curvature of the cover 
sheets does not conform with the curvature assumed by the ex- 
treme fibers of the spars when the beam is loaded. Formulas are 
derived for rib loads and for secondary sheet bending stresses as 
function of the primary sheet stress, the rib spacing, the beam 
height and the cover sheet thickness. With regard to the 
secondary stresses, the flexibility of ribs and spars, which in 
principle can be taken into consideration, is neglected so that the 
values obtained for the stresses are conservative. It is found that 
secondary sheet bending is largely restricted to narrow edge zones 
adjacent to ribs and spars. 


INTRODUCTION 


7 PURPOSE Of this article is to analyze the magni- 
tudes of a group of secondary stresses in the cover 
sheets of thin-walled box beams. The effects that will 
be considered in what follows can be described qualita- 
tively by saying that unless suitably stiffened the curva- 
ture of the cover sheets of box beams does not conform 
with the curvature assumed by the extreme fibers of the 
spars. The problem is analyzed here only for a cover 
sheet under tension, although a similar effect takes 
place in a cover sheet under compression before buck- 
ling sets in. The main results of the following develop- 
ments are analytical expressions for the cover sheet 
bending stresses and for the rib loads due to straighten- 
ing out of the cover sheets in between ribs and spars as 
functions of load intensity, rib spacing, beam height 
and cover sheet thickness. 

While in many practical cases the calculated second- 
ary stresses of this type will be small, it is shown by 
considering a numerical example that designs are pos- 
sible for which the magnitudes of the secondary stresses 
are an appreciable fraction of the magnitude of the pri- 
mary stress. . 

One of the conclusions from the analysis is the fact 
that the magnitude of the bending stresses obtained 
here may be such that experimental results concerning 
the effect of shear lag in cover sheets of box beams are 
appreciably distorted by the superposition of the 
straightening out effect. 


FORMULATION OF THE PROBLEM 


Consider a portion of a box beam between two adja- 
cent ribs. The rib distance should’ be small enough 
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Diagram of beam section in hetween two ribs: (a) unde- 
formed; (b) deformed. 


Fie. 1. 


that no appreciable change of bending moment takes 
place in the region considered (Fig. 1). 

According to elementary beam theory, the neutral 
axis of the beam and all fibers parallel to it become 
curved. The length R of the radius of curvature is 
given by 


1/R = M/EI (1) 


The spanwise stress in the cover sheets is given by 
a = (M/Dh (2) 


where h is the distance of the cover sheet plane from the 
neutral axis of the beam. 

On the tension side this stress tends to straighten out 
the sheet and may succeed in doing so except for the 
portions of the sheet adjacent to the spars. 

Since the straightened-out sheet fibers are shorter 
than the fibers attached to the slightly curved spars, 
a decrease toward the middle of the sheet of the direct 
spanwise stress ¢ will take place. It will be shown that 
in practical cases this decrease is exceedingly small and 
may safely be disregarded. 

On the other hand, it will also be shown that the 
regions of the sheet which are not completely straight- 
ened out are extremely narrow and thus assume an 




























SECONDARY STRESSES IN THIN-WALLED BOX BEAMS 


appreciable amount of secondary curvature, with the 
consequence that there may occur appreciable bending 
stresses in a narrow edge zone adjacent to ribs and 
spars. 


THE SECONDARY DEFLECTION OF THE COVER SHEET 
DvuE TO STRAIGHTENING-OUT AND THE RESULTANT 
REDUCTION IN Direct STRESS 


Introducing a system of coordinates x, y in the plane 
of the cover sheets, x spanwise and y transverse, the 
distance between curved and completely straightened 
sheet can be written as 


wo = nll — 4(x?/))) (3) 


where 7 represents the maximum secondary deflection 
and / the distance between ribs. 

Fig. 1b shows that 7 may be expressed in terms of h 
and R, 


h+R—1 = (h + R) cos (1/2R) (4) 


Since 


(l/2R)<1-°-h<«R 


the following formula can be written instead of Eq. (4) 
n = P/8R (5) 


Expressing R in terms of ¢ by means of Eqs. (1) and 
(2), there follows for the maximum secondary deflection 


n = (P?/8h)(o/E) (6) 


The reduction of the primary stress o is given by the 
compressive stress Ac, which reduces the length of the 
curved sheet to the length of the straightened-out sheet. 
The difference in length A/ of curved and straight sheet 
is given by 


Al = fu) V1 + diep/dx)? dx — 1 





(7) 


the integral in Eq. (7) being the analytical expression 

for the arc length of the curve w(x). Since (du/dx) 

< 1 Eq. (7) can be simplified to 

Al = (1/2) Ju), (dwe/dex)*dx (8) 

In terms of this change of length the compressive 
stress Ao is given by 


Ao = E(Al/l) (9) 


Introducing into Eq. (8) the function wp from Eq. (3) 
there follows first Al in terms of 7 and J, 


Al = (8/3)(n*/) (10) 
and hence Eq. (9) becomes 
Ao = E (8/3)(n°/?) (11) 


Introducing finally into Eq. (11) 7 from Eq. (6), there 
results for the relative reduction of stress due to straight- 
ening out of the sheet 


Aa/o = (1/24) (1/h)*(¢/E) (12) 
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Since o/E is not greater than '/0, it follows that for Ag 
to be more than '/j) of ¢ the distance between ribs has 
to be at least 50 times larger than the distance / be- 
tween cover sheet plane and neutral axis of the beam, 
which condition will certainly never even be approxi- 
mated. 

Thus, Eq. (12) shows that there is no need to take 
into account the deformation of the beam cross section 
for the calculation of the direct cover sheet stresses ex- 
cept in cases outside the range of practical interest. 


THE MAGNITUDE OF THE R1iB LOAD DUE TO 
STRAIGHTENING OUT OF THE SHEET 


Since the cover sheets of two adjacent bays, when 
straightened out, are inclined to each other, a com- 
ponent of the direct sheet stresses will be acting in the 
plane of the ribs. The magnitude of this component 
actually depends on the flexibility of the ribs. Here a 
conservative estimate of this magnitude is obtained on 
the basis of neglecting the flexibility of the ribs. 

From Fig. 1 it is seen that the angle of inclination.be- 
tween two adjacent straightened-out sheets is of magni- 
tude //R. The resultant cpg of the two inclined stresses 
o acts in the plane of the rib and is of magnitude 


or = (l/R)o (13) 

With R taken from Eqs. (1) and (2), og can be written 
in the form 

or = (l/h)(o?/E) 


To obtain the load ppg acting on a unit length of the rib, 
or has to be multiplied by the cover sheet thickness ?; 
hence, 


(14) 


be = ot(l/h)(o/E) (15) 


This load intensity pz will be acting over that portion 
of the width where complete straightening out takes 
place. It will be shown that this is the case all along 
the rib except in a narrow edge zone adjacent to the 
spars, so that it is only slightly more conservative to 
assume that the load intensity pg is uniform all along 
the rib. Eq. (15) indicates that the rib load intensity 
due to the straightening-out effect is directly propor- 
tional to the ratio of rib distance and cover sheet-neutral 
axis distance and proportional to the square of the di- 
rect sheet stress ¢. A numerical example will be con- 
sidered later. 


DETERMINATION OF THE SHEET BENDING STRESSES 


The determination of the bending stresses may be 
effected by considering the sheet as a thin plate having 
a slight initial deflection of magnitude w(x). The 
plate is acted upon by spanwise edge normal forces of 
magnitude ot per unit of edge length (Fig. 2). For this 
kind of loading the differential equation for the final 
deflection w, of the plate is known to be of the form 


V°V *w, = (ot/D)(0*w;/dx*) (16) 
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Fic. 2. Diagram for plate problem. 


where 

V? = (0?/0x") + (07/dy*) 
As boundary conditions are prescribed, along the ribs 

+1/2 w, = 0 Ow,/Ox = Ow/dx (18) 


and along the longitudinal supports, the distance be- 
tween which may be designated by 20, 


To make the boundary conditions of the problem homo- 


geneous, introduce the difference between the initial and 
final deflection 


D = Et8/12(1 — »)_ (17) 


x= 


(19) 


y= +h WwW, = WW 


W = Wy — Wy (20) 


into the system of Eqs. (16) to (19). There follows for 


w the following system of equations: 





eat (= i a) 21 
= D \Ox? ox? 1) 
x= +1//2 w=0 w/%x =0 (22) 
y= +b w=0 W/dy=0 (23) 


It should be noted that the boundary conditions of 
the problem are formulated as if the support offered to 
the sheet by ribs and spars were of an unyielding nature. 
This is certainly an idealizing assumption and the exist- 
ing flexibility of the supports will reduce the actual 
stresses below those calculated on the basis of rigid 
supports. There is then obtained a conservative esti- 
mate ouly of the magnitude of the effect under con- 
sideration, with the possibility left open to improve the 
accuracy of the results by modifying appropriately the 
conditions of support of the plate. 

The exact solution of the clamped-edge plate problem, 
Eqs. (21) to (23), offers considerable difficulty. 

It should be observed, however, that the sheet bend- 
ing along the longitudinal supports is dependent mainly 
on the conditions of support along the longitudinals and 
on the conditions along the ribs only to the extent that 
vanishing of the deflection there is essential, whereas the 
sheet bending along the ribs is due to the slope and de- 
flection restraint offered by the ribs, while the support 
offered by the longitudinals is restricted to longitudinal 

edge zones (as will be shown). Therefore, the following 
simplification can be introduced. The problem is sub- 
divided into two modified problems, one giving the 








1942 





stresses adjacent to the ribs (problem 1) and one giving 
the stresses adjacent to the spars (problem 2). 

In problem 1 the support offered by the spars is 
neglected, and with that the transverse variation of the 
deflection, so that the sheet is assumed to behave like 
an initially slightly curved beam acted upon by axial 
tension. 

In problem 2 the conditions of restraint at the ribs 
are assumed to be those of moment-free support. 


Determination of the Stresses Adjacent to the Ribs 


Considering problem 1, the differential equation and 
boundary conditions reduce to 


d‘*w ot d*w _ at 8n 
dxt Ddx®? DF 
+1/2 w=0 dw/dx =0 


(24) 








x= 25) 


where the value of w has been taken from Eq. (3). 
Putting as abbreviations 


i... 12(1 — v?) (26) 


(26a) 
the solution of Eqs. (24) and (25) becomes 


x? 
w=a[(1-42)—2 


Thus the final sheet elevation w; is given by the ex- 
pression 





cosh » — cosh kx 
— | (27) 
‘p sinh pw 


cosh nw — cosh kx 





MW=—-Hy—-wWw= 2n (28) 


uw sinh pu 
For the secondary bending stress in the extreme sheet 
fibers, which is due to the relative deflection w, there 
follows from Eq. (27) 


E tdw 


:, => = aW“— = + 


1 — pv? 2 dx? 








E t (- 8 2k? cosh =) 
i— 93" }? p sinh p 
which may be written in the form 

E Atn (" cosh kx 


1—- - 





_ 1) (29) 


Pei sinh p 
Inspection of Eq. (29) reveals that o, decreases rapidly 
from the edges x = +//2 when u > 1, which is also 
found to be the condition for occurrence of appreciable 


bending stress. Writing then 
” “(1-75 
cosh kx/sinh uw ~ e W/2 (30) 
there follows from Eq. (29) 
E 4tn aig tone. 
o; = “I-97 * u( 1/2 (31) 
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Eq. (31) shows that the maximum stress occurs di- 
rectly adjacent to the ribs and decreases rapidly so that 
at a distance 


6, ~ 1.2(1/u) (32) 


it amounts to only 10 per cent of the edge stress. In- 
troducing the value of u from Eq. (26) and the value of 
» from Eq. (6), there follows for the maximum bending 
stress adjacent to the ribs and for the width of the edge 


zone 
o, 0.75 1 [eo 
Sm wg Ot ae (38) 
E 
5, ~ 0.8¢ \. (34) 
Co 


Eq. (33) indicates that o, increases more rapidly than 
o and that it is directly proportional to ratio of rib dis- 
tance / and distance / from sheet to neutral axis, which 
ratio occurred previously in Eq. (15) for the rib loads 
pr. Eq. (24) shows that the width of the edge zone is 
proportional to the sheet thickness and decreases with 
increasing primary stress oc. 

It may be said that the value of the rib load intensity 
pr as given by Eq. (15) can also be obtained as the dif- 
ference in vertical shear between two adjacent ribs— 
that is, from the expression 


Vopr = —D(@*w/dx*),=1/2 


Determination of the Stresses Adjacent to the Spars 

To obtain values of the bending stresses along the 
spars, problem 2 is to be solved, which consists of the 
differential equation, Eq. (21), the boundary conditions, 
Eqs. (23), and the modified boundary conditions 


x= +//2 w=0 M,=0 (35) 


A suitable solution of the differential equation is ob- 
tained by assuming 


w = cos ax fy) +0 (1 - 42) a=n/l (36) 


This satisfies the displacement condition in Eqs. (35) 
exactly and the moment condition very nearly, so that 
these conditions are considered to be taken care of. 

Introducing Eq. (36) into the differential Eq. (21), 
there follows for the function f(y), 


d‘f , of ( oe ) 
—_ — — = =0 37 
dy! wal ° ta ate f (37) 
which with the notations of Eqs. (26) may be written 
d'f Uf ( 4 ‘) 
—_— — — = =0 38 
Zi stellt oe )s (38) 
Since, as before, 1 > 1 within the range of interest, the 


solution of Eq. (38) is very nearly the same as the solu- 
tion of the following simpler equation, 


(d*f/dy*) + (4/x*)watf = 0 (39) 
Putting 
aty?/x? = 4 (40) 
the solution of Eq. (39) is of the form 
f = e~7- (C sin vy + Cocos vy) + 
e’©-» (Cs sin yy + Cy cos yy) (41) 
This introduced into Eq. (36) represents the solution for 
w in which the constants of integration C,, are to be de- 
termined by means of the boundary conditions Eqs. 
(23). 
The form of the function f(y) indicates that its effect 
will be restricted to zones in the immediate neighbor- 


hood of the edges y = +3, and it is possible to simplify 
the solution by considering separately the effects along 


the two edges. Considering the edge y = 6, C3 = Cy = 
0 may be put so that 
w= cos ax e~?°’-”) (C, sin vy + Crcos yy) + 
n [1 — 4(x°/1?)] (42) 
The boundary conditions, Eqs. (23), are sufficiently 
nearly satisfied in view of the similarity of the curves 
represented by cos ax and by [1 — 4(x?//*)] by pre- 
scribing w(0, b) = 0 instead of w(x, 5) = 0 so that the 
equations for the constants of integration become 
w(0,b) = 0 C,sin yb + C. cos yb = —7n (43) 
(Ow/Oy)y-1 = 0 C,(sin yb + cos yb) + 
C:(cos yb — sin yb) = 0 (44) 
Solving for C; and C; and substituting in Eq. (42), 


w=n {cos axe~7-” [cos y(b — y) + 


sin y(b — y)] — (1 —4 =) (45) 


For the bending stress normal to the spars in the 
extreme sheet fibers there follows, in accordance with 


plate theory, 
E : (= ove) E t dw 


+ —— — | — y— 
1 — v® 2 \dy’? Ox? 





Oy = 


the term proportional to v being negligible because ¥ is 
large compared with a. Hence, 
E ft 


T—y2"* 





1= F 


{22 cos ax e~7°- [cos y(b — y) — sin y(b — y)} } 
Taking of Eq. (40) which has the value 


re ayen YEE 
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and 7 from Eq. (6) so that 


2 


t 2=> ft— 
nY Sh 





there follows, if one also puts 
V3 
8V 1 — rv? 


= 0.71 VW 3r2(1 — v?) = 2.28 


2 A Ber — sin |: e =o ) 
es (2. (2) (2 2) (49) 


The stress o, takes on its maximum value o, adjacent to 
the spars and halfway in between the ribs. From Eq. 
(49) it is seen that 


a, = 0.7lo(l/h)Vo/E 


Comparison of o, with o, of Eq. (33) shows identical 
analytic expressions except for the value of the numeri- 
cal coefficient, which in o, is about 20 per cent larger 
than in o;. 

Eq. (49) indicates that the o, curve has the character 
of a damped oscillation, however, so rapidly damped 
that the values of o, beyond the first half wave are small 
fractions of as. 

The width 6, of the edge zone can here be defined as 
the distance from the spar to where a, becomes zero for 
the first time. From Eq. (49) it follows that this is the 
case when 


(50) 


2.28Wo/E (b — y)/Vit = 2/4 
whence 
5, = 0.35Vil WE/o (51) 
Comparison of Eq. (51) with Eq. (34) shows that there 
are different analytic laws for 6, and for 6, but that in 
both cases the occurrence of the sheet thickness ¢ as a 
factor indicates the narrowness of the edge zone. It is 
because of this narrowness that the secondary sheet 
curvatures and, with that, the sheet bending stresses 
in the edge zone are rather larger than would otherwise 
be the case. 

With the help of the solution obtained it is also pos- 
sible to determine the forces between sheet and spars 
due to the straightening out of the sheet. These forces 
are the vertical reactions corresponding to the plate 
deflection as given by Eq. (45). According to plate 
theory the load intensity acting on the spars is here 

o*w o'w 
b, = —D ( 5+ @- 5 ) (52) 
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Because y is large compared with a the second term on 
the right side of Eq. (52) may be neglected. Introduc- 
ing w from Eq. (45) into the simplified Eq. (52) there 
follows 

(53) 


pb; = —DA4y'n cos ax 


and introducing the values of 7, D, and y from Egs. 
(6), (17), and (47) there results, finally, 


vi (z)" x 
—0.54c0t — | — 
0.540 a \E COS © 7 


= (54) 
Comparing this result for the spar load intensity p, with 
Eq. (15) for the rib load intensity pp, it is seen that, in 
contrast with the expressions for the sheet bending 
stress, p, and pz are given by different analytic expres- 
sions. The important difference between p, and pp is 
that a factor / in Eq. (15) is replaced by a factor ~/#/ 
in Eq. (54) so that p, is of a smaller order of magnitude 
than pp. 


A NUMERICAL EXAMPLE 


To obtain an impression of the magnitude of the 
secondary stresses, a beam section of the following dimen- 
sions may be taken: rib distance / = 10 in.; distance 
of cover sheet from neutral axis h = 5in.; sheet thick- 
ness ¢ = 0.1 in.; cover sheet stress ¢ = 105 Ibs. in.-?; 
modulus of elasticity E = 107 Ibs. in.—? 

There follows: ; 

(1) From Eq. (15) for the rib load per inch, 


Pre = 200 Ibs. in.-! 


(2) From Eq. (33) for the secondary sheet bending 
stress adjacent to the rib, 


o, = 0.180 = 18,000 Ibs. in.-? 


(3) From Eq. (48) for the secondary sheet bending 
stress adjacent to the spar, 


o, = 0.14¢ = 14,000 Ibs. in.~? 


(4) From Eq. (34) for the width of the edge zone at 
the rib, 


5, = 0.8t = 0.8 in. 


(5) From Eq. (49) for the width of the edge zone at 
the spar, 


6, = 1.17 10¢ = 1.1 in. 
(6) From Eq. (54) for the spar load intensity, 
Pb. = —35 cos (rx/I) Ibs. in.-! 


Finally, one may calculate the value of the maximum 
initial deflection 7 from Eq. (6), 
n = 0.025 in. = 0.25¢ 
This result for » confirms the supposition that it is 


appropriate to apply the small deflection theory of plates 
to the problem at hand. 











SUMMARY 


Theoretic equations are developed to determine what concen- 
trations of oxygen must be breathed at any altitude to maintain 
the “oxygen equivalent” of air at given altitudes—i.e., the same 
alveolar oxygen pressure. The partial pressures of water vapor 
and carbon dioxide are taken into account, as well as the varia- 
tion of the pressure of carbon dioxide with the alveolar oxygen 
pressure. 

Six curves are given on a chart showing the oxygen concen- 
trations required at any altitude to maintain the oxygen equiva- 
lent of air at sea level and at altitudes of 5,000, 10,000, 15,000, 
20,000, and 25,000 ft. Altitudes of 33,800, 40,000, and 46,000 
ft. with pure oxygen are equivalent to air at sea level, 10,000 and 
20,000 ft., respectively. 

It is shown that these equivalents are nearly independent of 
the assumptions as to the state of physiologic adaptation and 
should be applicable to all individuals. 





[ IS GENERALLY RECOGNIZED that most of the physio- 

logic dangers of high altitude are due to lack of 
oxygen and may be largely avoided by the use of 
supplemental oxygen to maintain normal oxygen 
pressure in the lungs. The main exception is the 
danger from aeroembolism, which arises from gases 
coming out of solution in the body. Both qualitative 
and quantitative studies have been made on aero- 
embolism, ? and a great deal of experimentation has 
been done on the qualitative aspects of anoxia." ® 
However, data on the quantitative aspects of preven- 
tion of anoxia are still rather scarce. 

The important question—What percentages of 
oxygen are required in the inspired air at various 
altitudes to maintain an alveolar oxygen pressure 
equivalent to that obtained with a given oxygen per- 
centage at a given altitude?—has not been experi- 
mentally investigated in detail. It is the purpose of 
this note to present the derivation of a formula and a 
chart based on it for theoretic determination of such 
respiratory equivalents. 


RESPIRATORY EQUATIONS 


The respiratory equations must be deduced from 
considerations of the effect of water vapor and CO, on 
the partial pressure of oxygen in the inspired air. 

For convenience, the symbols used in the develop- 
ment will be defined here. 
ventilation, or volume rate of breathing, including 


atmospheric air and added oxygen, if any, in 
liters per minute, at ambient pressure and 0°C. 


V = 





B = ambient atmospheric pressure, mm. Hg 

POA = partial pressure of oxygen in dry inspired air, 
mm. Hg 
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POL woz. = Maximum oxygen pressure in upper part of lungs, 
mm. Hg 

POL = partial pressure of oxygen in the lung alveoli, 
mm. Hg 

PCL = partial pressure of CO, in the lung alveoli, mm. Hg 

c = POA/B = fractional concentration of oxygen 

PO = decrease in oxygen pressure in lungs due to 


absorption, mm. Hg 


The maximum oxygen pressure in air moistened 
by passage through the respiratory tract is 


POL mar. = POA[(B — 47)/B] 
At sea level 
POL nar. = [(21/100) X 760] X [(760 — 47)/760] = 150 


It may be noted that this is independent of the hu- 
midity of the inspired air, moisture being supplied by 
evaporation as required so that the air is finally satu- 
rated at nearly body temperature. 

The effect of CO. must also be considered explicitly. 
The inspired air is diluted with water vapor, of pressure 
47 mm. Hg, largely while being drawn through the 
nasal passages; but, before reaching the absorbing 
tissues of the alveoli, it is further diluted by COzs, 
of a pressure normally about 40 mm. Hg. Thus, the 
maximum oxygen pressure before absorption would be 


POL = POA[(B — 47 — PCL)/B] (1) 


which, at sea level, equals 140 mm. Hg. 

If only oxygen, CO, and water vapor are present, 
POL in Eq. (1) is also the effective oxygen pressure at 
the absorbing alveolar walls. 

If nitrogen or other diluent is present, local absorp- 
tion of oxygen will effectively reduce the partial pres- 
sure of oxygen. Analysis of the last part of the air 
in a forced exhalation (alveolar air) indicates that the 
oxygen pressure in the alveoli averages about 100 mm. 
Hg for normal men breathing air at sea level (Fig. 1). 
This difference would indicate a pressure reduction 
due to absorption of about 40 mm. Hg; but, for the 
moment, call this decrease merely APO. 

To generalize these equations for all altitudes, for 
all oxygen concentrations, and for all persons, it is 
necessary to use the values for alveolar CO: pressure 
and of the oxygen consumption, characteristic of the 
person, and allow for their possible variation. 

The equation 


POL = POA([(B — 47 — PCL)/B] — APO (2) 


is perfectly general but useless until values for PCL 
and APO are given. PCL may vary from 20 to 40 
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Alveolar gas pressures for persons at various oxygen 
pressures or altitudes. Curves (1) and (4) are for men acclima- 
tized at the indicated altitudes; (2) and (5) are for men taken 
from sea level to the indicated altitudes; (3) and (6) are for per- 
sons breathing air-oxygen mixtures to maintain normal O,2 and 
CO, pressures in the lungs. 


Fic. 1. 


mim. Hg, depending on the altitude to which the person 
is accustomed. Since the effects of high altitudes 
on persons who live near sea level are of the most 
concern, PCL may usually be taken as equal to 40. 
However, one effect of anoxemia is to increase the 
breathing volume, which washes the CO: out more 
rapidly and lowers its partial pressure, so that under 
conditions of plain anoxia, the approximation 


PCL = 40(Vo/V) (3) 


might be used, where Vo is the normal, and V is the 
increased, rate of ventilation. At sea level APO is 
about 40 mm. Hg when air is breathed and, as has 
been mentioned above, will be zero when pure oxygen 
is breathed. As a possible law of variation with 
concentration, let 


APO = 50[1 — (POA/B)] (4) 


which satisfies the two conditions given with a linear 
dependence on concentration. The value of 50 in this 
equation is only an average figure; it depends on the 
respiratory quotient and on the ratio of ventilation 
to oxygen consumption characteristic of the individual. 

As mentioned earlier, the presence of CO, in the 
inspired air stimulates breathing, as does anoxia also 
to a lesser extent. With increased breathing, the 
ratio of oxygen absorbed to the volume breathed may 
decrease. The effect of such increased breathing may 
therefore be taken into account by assuming 


APO = 50[1 — (POA/B)}(Vo/V) (5) 
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which gives, when inserted into Eq. (2), 


POL = POA[(B — 47 — PCL)/B] — 
50[1 — (POA/B)](Vo/V) (6) 


In this general equation, PCL will have nearly its 
normal value of about 40 mm. Hg when any reasonable 
amount of CO, is present in the inspired air and will 
otherwise have values as given by Eq. (3). 

In aviation one is generally concerned with persons 
who fly from near sea level to high altitudes. The 
physiologic reactions of such persons are different from 
those of persons acclimatized to higher altitudes and 
will depend also on whether oxygen, air, or oxygen-air 
mixtures are breathed. 

Fig. 1 shows the average values of alveolar CO: and 
O» pressure in men acclimatized at various pressures or 
altitudes (lines 1 and 4, adapted from Haldane‘) and 
for men taken rapidly from sea level to the indicated 
pressures or altitudes, breathing air only (lines 2 and 5 
adapted from McFarland* and approximated as straight 
lines). The horizontal lines 3 and 6 indicate the 
normal sea-level values, which may be maintained up 
to some limiting altitude, by breathing the proper 
mixture of supplemental oxygen and air. What this 
limiting altitude is and what mixtures are required will 
be developed below. 

The lines 2 and 5 in Fig. 1 represent transient condi- 
tions, in one sense; after some time at a given altitude, 
the representative points approach the curves 1 and 4 
as acclimatization occurs, accompanied by further 
increase in the normal ventilation rate. 

The rapidity with which adaptation occurs varies 
widely among individuals; a few hours are sufficient 
for considerable accommodation in some, while days 
or even weeks are required for others. Because of 

individual variations it is impossible to derive an 
equation based on averages that will satisfactorily 
predict the physiologic state of all individuals under 
anoxic conditions. However, a great deal may be 
accomplished with considerable generality, if the aim 
is changed only slightly. Let it be attempted to 
derive equations relating oxygen concentrations to 
altitude, such that one can determine what concen- 
tration of oxygen is required at a given altitude (or 
ambient pressure) to maintain the respiratory equiva- 
lent of a given concentration at some other altitude. 

Eq. (6) has too many variables as it stands; but, 
if some relation is assumed between PCL and POL or 
PCL and B, an equation of the desired form is ob- 
tained. The corresponding O, and CO, curves in 
Fig. 1 supply three different relations between PCL 
and POL: 


(a) from line 6: No acclimatization, 
PCL = 40 
(b) from lines 1 and 4: Partial acclimatization, 


PCL = 34 + 0.6POL (8) 


(c) from lines 2 and 5: Complete acclimatization, 


PCL = 12 + 0.28POL 


(9) 


















There is no evidence that anyone fails to adapt to 
anoxia to at least some degree, so that (a) is clearly an 
extreme; average individuals may be expected to 
approximate (b) or (c), depending on the time of 
exposure to low oxygen concentrations and individual 
adaptability. 

It is probable that for short-time exposure to various 
altitudes the partial pressure of CO, in the lungs does 
not vary linearly with the ambient pressures; some 
curve intermediate between the lines shown in Fig. 1 
may be more representative for most individuals. 
However, the straight lines are easier for computation, 
and the results based on them will serve as limits. 

Substituting these three relations in Eq. (6) and 
putting POA/B = C, and from Eq. (3), Vo/V = 


PCL/40, the following three equations result. For 
the different degrees of adaptation: 
C = (POL + 50)/(B — 37) (10a) 


C = (1.075POL + 42.5)/(B — 38.5 + 0.015POL) (10b) 
C = (1.35POL + 15)/(B — 44+ 0.07POL)  (10c) 


If Cis put equal to 0.21 and the value of B correspond- 
ing to various altitudes is inserted, the values of POL 
in air at selected altitudes can be obtained on the three 
assumptions. Using the values of POL thus deter- 
mined, relations between C and B are obtained, giving 
the concentrations necessary at any altitude to maintain 
the POL values corresponding to the selected altitudes 
in air. 

For comparison of the extremes, the alveolar oxygen 
pressures and the highest altitudes at which they can 
be maintained (C = 1) are listed in Table 1, computed 
on the basis of Eqs. (10a), (10b), and (10c). 


TABLE 1 
Alveolar Oxygen Pressures at Various Altitudes, and Highest 
Altitudes at Which They May Be Maintained with 100 Per Cent 
Oxygen, Computed on the Bases of Eqs. (10a), (10b), and (10c) 


OXYGEN REQUIREMENTS AT HIGH ALTITUDES 








Equivalent 
Altitude in 
Limiting Altitude with 


Air (21% Alveolar Oxygen 

Oxygen) Pressure (mm. Hg) 100% Oxygen (1,000 Ft.) 

1,000 Ft. (10a) (10b) (10c) (10a) (10b) (10c) 
0 102 102 102 33.8 33.8 33.8 
5 75 if 81 37.0 37.0 37.0 
10 52 55 64 40.2 40.1 39.9 
15 30 37 49 43.5 43.3 43.0 
20 16 21 37 46.6 46.4 45.9 
25 1.5 8 26 49.7 49.5 48.8 





In spite of the great differences in the values of POL 
computed on these divergent assumptions as to degree 
of adaptation of the individual, the limiting altitudes 
are only slightly affected. It may be concluded that 


any one of the equations given furnishes a good repre- 
sentation of respiratory equivalents—that is to say, 
for instance, that whatever may be the respiratory 
state of any individual breathing air at 20,000 ft. 
breathing pure oxygen will maintain this state at about 
46,000 ft. 
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To compute circulatory conditions and the general 
physiologic state from respiratory factors requires 
specific assumptions as to the mode of variation of 
POL, but this is not the object here. 

One point should be mentioned to guard against 
confusion. The limiting altitudes in Table 1 are seen 
to be slightly higher for the case where no acclimatiza- 
tion is assumed. It must be remembered that these 
limiting altitudes are those for which the respiratory 
condition is the same as for the specified altitude, when 
breathing air, for the degree of acclimatization as- 
sumed. Thus, while the respiratory state of a poorly 
adapted individual is the same at 43,300 ft. with 100 
per cent oxygen as at 15,000 ft. with air, it is clearly not 
so physiologically favorable as that of an individual 
acclimatized at 15,000 ft. The unacclimatized person 
may be severely handicapped, while the perfectly 
adapted person can carry on normal physical activities 
under either condition. Even though the latter can 
maintain his higher POL completely only up to 43,000 
ft. according to Table 1, the decrease in his POL 
in 300 ft. would be only a few millimeters of mer- 


cury. 
CHART OF OXYGEN EQUIVALENTS 


Since the condition assumed in Eq. (10c) is the one 
toward which average persons tend with time, and 
particularly since it gives the more conservative 
limiting altitudes, it has been used in computing the 
chart of Fig. 2, which shows the concentration of 
oxygen required at any altitude to maintain the re- 
spiratory equivalent of air at sea level and at altitudes 
of 5,000, 10,000, 15,000, 20,000, and 25,000 ft. 

100 


OXYGEN REQUIREMENT CHART 
NATIONAL BUREAU OF STANDARDS 
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Fic. 2. Oxygen Requirement Chart. Curves show the per- 
centage of oxygen necessary at other altitudes to maintain alveo- 
lar O2 pressures equivalent to breathing air at sea level, 5,000, 
10,000, 15,000, 20,000, and 25,000 ft. 


The sea-level curve—the one most generally used— 
is not affected at all by the differing assumptions as 
to variation of POL and PCL in anoxia, while the other 
curves would spread out only slightly more on the basis 
of the other equations, as indicated in Table 1. 
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For comparison with the other extremes, a broken 
curve based on Eq. (10a) is shown for the 25,000-ft. 
equivalent, which gives the greatest spread. The 
fact that the difference between the extremes is slight 
is a justification for the approximations made in mathe- 
matical representation of the extremes. 

No alveolar oxygen pressures are shown on the 
chart—it may be used to determine the concentrations 
required to maintain oxygen equivalents regardless 
of the exact value of POL which any particular indi- 
vidual might have under an assumed condition. As 
pointed out above, the concentrations in the inspired 
air required to maintain the oxygen status quo in the 
lungs at other altitudes depend only slightly on the 
exact value of POL. Any computation of POL 
values for aviators should be based on Eq. (10b); 
however, since it gives the more conservative values 
for this quantity. 

From this it must not be concluded that the total 
oxygen requirements are independent of the POL 
value for the individual—far from it. The increase 
in POL gained by acclimatization or adaptation is 
accompanied by a considerable increase in ventilation 
so that the total volume of oxygen breathed is corre- 
spondingly increased. If the oxygen in the exhaled 
breath is wasted, the total oxygen to be supplied is 
therefore also increased. 

On the auxiliary ordinate scale of the chart are shown 
the percentages of added oxygen in oxygen-air mixtures 
required to give the desired oxygen concentrations. 
The total oxygen fraction, C, in the inspired air is 
related to the oxygen added, expressed as a fraction of 
the total volume, C,, by the simple equation 


C = 0.79C, + 0.21 (11) 


This chart is useful in determining the sufficiency of 
various oxygen concentrations at various altitudes and 
in estimating the safety factor in oxygen supply. A 
similar chart based on Eq. (10b) has been used for 
several years at the National Bureau of Standards and 
the Bureau of Aeronautics, Navy Department, in 
specifying and evaluating the performance of oxygen 
supply systems. 

It is seen that the percentage of oxygen required to 
maintain a given respiratory equivalent increases more 
rapidly at the higher altitude. A 50-50 mixture of 
added oxygen and air will give a sea-level respiratory 
equivalent at nearly 25,000 ft., while 100 per cent O 
is required for this equivalent at 33,800 ft. Also, the 
effective respiratory altitude increases much more 
rapidly than the actual altitude when breathing oxygen 
or oxygen-air mixtures. Ascent from 33,800 to 
46,300 ft. (12,500 ft.) with 100 per cent oxygen is 
equivalent to a change from sea level to 20,000 ft., 
breathing air. 

As another example of the use of the chart, suppose 
that in flight at 25,000 ft. breathing conditions equiva- 
lent to 10,000 ft. in air are desired. Following the 


third curve—for air at 10,000 ft.—up to the inter- 
section with the 25,000-ft. ordinate, the ordinate scale 
shows that 43 per cent oxygen is needed, which, by 
the auxiliary scale, can be obtained by an air-oxygen 
mixture with 28 per cent added oxygen. 

Since normal individuals do not lose noticeably in 
efficiency up to 10,000 ft. in air, with pure oxygen 
they should be equally able at 40,000 ft. (upper end of 
10,000-ft. curve), as far as the effects of oxygen pressure 
itself are concerned. 

To be sure, aeroembolism may occur at this altitude 
and even cause complete disability, although this is 
quite apart from the oxygen pressure in the lungs. 
Also, there may be other effects of reduced pressure 
which affect the physiologic state, so that one may 
not immediately identify “‘oxygen equivalents’ with 
circulatory equivalents, or general fitness equivalents, 
until or unless experimental data become available to 
warrant it. 

In recent experiments by Dill and Hall* on three 
men breathing pure oxygen at 44,000 ft., arterial 
saturations equivalent to about 18,000 ft. in air were 
found. The chart indicates that air at 18,000 ft. is 
the ‘‘oxygen equivalent” of 100 per cent oxygen at an 
altitude of about 44,800 ft. In the light of the fore- 
going remarks and the small number of experimental 
data, no attempt will be made to assess here the signifi- 
cance of this difference. 


EFFECT OF WATER VAPOR AND CO; 


In the development of the equations on which the 
chart is based, inspired mixtures of oxygen and dry 
air only were implied. This restriction is clearly 
unnecessary. Since the oxygen is diluted by water 
vapor and carbon dioxide before reaching the zone of 
absorption, it makes no difference whether the dilution 
occurs before or after inspiration. Thus, the addition 
of water vapor to the dry oxygen-air mixtures indicated 
by the chart will not change the respiratory equiva- 
lent. Likewise, the addition of small amounts of 
CO, (partial pressure up to about 20 mm. Hg) to the 
mixture will not change the respiratory equivalent if 
the same partial pressure is added at all altitudes. 

The effect of added CO, is mainly to increase the 
ventilation. When breathing air or air-oxygen mix- 
tures this may result in better arterial saturation, as 
noted by Haldane! and by Dill. As indicated by 
Eqs. (5) and (6), increased ventilation will reduce the 
loss of oxygen pressure due to absorption (APO), 
thereby increasing POL. Also, the better ventilation 
may result in more thorough gas exchange in parts 
of the lungs which are otherwise poorly ventilated. 
However, in flight at the highest altitudes, anoxia 
alone causes greatly increased ventilation, APOL is 
already zero if pure oxygen is breathed, and the partial 
pressure CO, in the lungs is lowered. Under such 
conditions, breathing in enough CO, to raise the partial 
pressure of CO, in the lungs could only result in lower- 
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ing the effective oxygen pressure. On this basis CO: 
would be of doubtful value. This deduction is in 
general agreement with the experimental findings of 
Dill.* 


ADVANTAGES OF ACCLIMATIZATION 


As an example of the advantage that a person al- 
ready acclimatized at a high altitude would have in 
extreme high-altitude flight with oxygen, consider the 
state of a person acclimatized at sea level in comparison 
with that of a person acclimatized at 20,000 ft., when 
both are supplied with oxygen and taken to, say, 
49,000 ft. From Table 1, columns 2 and 3, it is seen 
that at 20,000 ft. the person from sea level will have 
an alveolar oxygen pressure of only 21 mm. Hg, as 
against 37 mm. Hg for the acclimatized person. 

But collapse occurs at just about this lower value of 
POL; so with 100 per cent oxygen the person from 
sea level would still be on the verge of collapse at 
46,000 ft., while the acclimatized person would be in 
his “normal’’ condition. To go much higher would 
be fatal for the one subject; the other might ascend 
another 3,000 ft.—or conceivably even 4,000 ft.— 
without losing consciousness. According to the chart, 
the air equivalent of 49,000 ft. with oxygen is 25,000 
ft., which is tolerable for persons acclimatized at 
20,000 ft. 

It is of interest in this connection that the record 
altitude attained in flight without a pressure suit 
(Donati, Italy, 1934) was to 47,358 ft. and was ap- 
parently very near the limits of the physical endurance 
of the pilot. It is not recorded whether he had the 
benefit of any prior acclimatization. 

Another theoretic speculation of some interest 
concerns the oxygen problem in attempted ascents of 


Mount Everest (elevation 29,000 ft.). On two of six 
expeditions altitudes above 28,000 ft. were achieved 
with and without the use of oxygen. No data on al- 
veolar air composition was obtained at this altitude, 
but a sample taken after several days at 22,700 ft. 
indicated acclimatization, with alveolar oxygen pres- 
sure of 38.8 mm. Hg, and CO: pressure of 19.3 mm. Hg, 
in fair agreement with the curves of Fig. 1. , 
The chart indicates that such climbers, already ac- 
climatized to 20,000 ft. or thereabouts, could maintain 
their accustomed respiratory conditions up to the 
summit by breathing 20 per cent added oxygen, whereas 
a person acclimatized at sea level would require about 
70 per cent, to maintain sea-level equivalent. How- 
ever, the breathing volume for the person accustomed 
to 20,000 ft. would be double the value normal for sea 
level. If the exhaled oxygen were wasted, the ad- 
vantages of acclimatization in reducing the required 
oxygen concentration would therefore be somewhat 
offset. Such considerations are of importance in 
determining what type of apparatus is best suited to 
particular uses. ; 
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rectly applicable to airplanes for private use. It is be- 
lieved that such a plane in the rather near future will 
present the following characteristics: (1) it will have 
outstanding vision, (2) it will neither spin nor wind up 
in a spiral, and (3) it will be stable and consistent in 
control throughout its entire maneuvering range. 


The controls will present no change from present con- 
ventional practice that would confuse a normal pilot 
with conventional training. Yet they will appeal to the 
average intelligent person who drives a car, who finds it 
easy enough to handle one more control, but who is 
confused by any further additions. 


The control wheel, looking like the wheel in a car, 
will steer in similar fashion—whether in the air or on 
the ground. Its fore-and-aft movement will coordinate 
all functions of lift. A pedal will provide drag or brake 
force either in the air or on the ground. The throttle 
will be preferably a foot type but with means of holding 
any desired amount of cruising power. For most pur- 
poses dual controls will be as superfluous as in an auto- 
mobile. 

At the same time, performance can be stepped up to 


an extent that will leave no doubt as to the superiority 
of airplane over automobile for cross-country trips. 
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It seems that Newton felt that the real solution of the 
resistance problem required finding the pressure dis- 
tribution in the entire fluid space. This he was not 
able to do since he did not possess the methods needed 
to go deeper into this problem. He felt that the pres- 
sure propagation decreased the resistance—decreased 
it, in fact, to nil in the case of a nonviscous fluid, as it 
was found later by d’Alembert. He could not see that 
the pressure propagation increased tremendously the 
normal force. He did not consider this case at all; 
apparently he did not connect the idea of flight with the 
concept of air resistance. 


There are two limiting cases, however, in which New- 
ton'’s law yields approximately correct results, because 
the physical picture is similar to that visualized by him. 
It was shown by recent theoretic and experimental in- 
vestigations that the lift coefficients of a plane airfoil of 


very small aspect ratio—i.e., of an inclined strip, whose 
dimension in the flight direction is large in comparison 
with its dimension normal to the flight direction—is 
given with fair approximation by Newton’s Eq. (2). 
In this case, the deflected air mass is approximately 
equal to the air mass that hits the airfoil.! 

The second case refers to the theory of drag for shells 
of extremely large Mach numbers (approaching infin- 
ity). In this case, the shock wave is close to the surface 
of the shell. Therefore, the pressure practically does 
not propagate into space and the effect of the wind 
stream is similar to that of an inelastic impact as as- 
sumed in Newton’s theory. 
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cases constitute the only source of information used in compiling 
scientific reference indexes. Abstracts printed in other jour- 
nals, especially foreign, in most cases, consist of summaries 
from printed papers. The summary should explain as adequately 
as possible the major conclusions to a nonspecialist in the subject. 
The summary should contain from 100 to 300 words, depending 
on the length of the paper. 


Sus-HEapINGs: Sub-headings should be inserted by the au- 
thor at frequent intervals. The work of editorial preparation 
will be simplified by the author providing many sub-headings. 
Owing to the breaking of columns and the insertion of illustra- 
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